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1 Introduction and Overview^ 



Previous studies [||] of the physics potential of LEP2 indicated that with the design luminosity 
of 500 pb~ one may get a direct measurement of the W mass with a precision in the range 
30 — 50 MeV. This report presents an updated evaluation of the estimated error on Mw based 
on recent simulation work and improved theoretical input. The most efficient experimental 
methods which will be used are also described. 

1.1 Machine parameters 

The LEP2 machine parameters are by now largely determined. Collider energy values and time- 
scales for the various runs, expected luminosities and errors on the beam energy and luminosity 
are discussed and summarized elsewhere in this report 0, Here we note that (i) collider 
energies in the range 161 — 192 GeV will be available, and (ii) the total luminosity is expected 
to be approximately 500 pb~^ per experiment. It is likely that the bulk of the luminosity will 
be delivered at high energy {y/s ^ 175 GeV). The beam energy will be known to within an 
uncertainty of 12 MeV, and the luminosity is expected to be measured with a precision better 
than 1%. 

1.2 Present status of Mw measurements 

Precise measurements of the masses of the heavy gauge W and Z bosons are of fundamental 
physical importance. The current precision from direct measurements is AMz = ± 2.2 MeV and 
AMvK = ± 160 MeV So far, Mw has been measured at the CERN |^ and Fermilab Tevatron 
[H) 0) m PP colliders. The present measurements are summarized in Fig. |I[ In calculating the 
world average, a common systematic error of ±85 MeV arising from uncertainties in the parton 
distributions functions is taken into account. The current world average value is 

Mw = 80.26 ± 0.16 GeV. (1) 

An indirect determination of Mw from a global Standard Model (SM) fit to electroweak data 
from LEPl and SLC ||^ gives the more accurate value 

Mw = 80.359 ± 0.051 toml GeV . (2) 

In Fig. m this range is indicated by dashed vertical lines. Note that the central value in (H) 
corresponds to Mh = 300 GeV and the second error indicates the change in Mw when Mh is 
varied between 60 GeV and 1000 GeV - increasing Mh decreases Mw. The direct measurement 
of Mw becomes particularly interesting if its error can be made comparable to, or smaller than, 
the error of the indirect measurement, i.e. < 0(50) MeV. In particular, a precise value of Mw 

^prepared by F. Jegerlehner, Z. Kunszt, G.-J. van Oldenborgh, P.B. Renton, T. Riemann, W.J. Stirling 



4 



indirect (1a limits) — > 




< — 


CDF 90 Q 






79.91 ± 0.39 






UA2 92 






80.35 ± 0.37 


o- 




DO (prelim.) 94 

o 

79.86 ± 0.35 






CDF 95 






80.41 ± 0.18 






Average 8o.26±o.i6 

— • 






85 MeV common error 
. . . . 1 , . , . 1 . , . , 1 , . , . 1 . . . . 1 . . . . 1 . , , , 1 




. 1 .... 1 



78.5 78.75 79 79.25 79.5 79.75 80 80.25 80.5 80.75 



Mw (GeV) 

Figure 1: Direct measurements of M-\j^ together with their average. The dashed lines correspond 
to the indirect limits (± 1 a) from a global fit in the Standard Model. 



obtained from direct measurement could contradict the value determined indirectly from the 
global fit, thus indicating a breakdown of the Standard Model. An improvement in the precision 
of the M-w measurement can be used to further constrain the allowed values of the Higgs boson 
mass in the Standard Model, or the parameter space of the Minimal Supersymmetric Standard 
Model (MSSM) . 

Standard Model fits to electroweak data determine values for Mw (or Mh), Mt, and 



as{Mz). The direct determinations of the top quark mass |^ give an average value of 
Mt = 180 ± 12 GeV. Fig. ^ compares the direct determinations of Mt and M-w with the in- 
direct determinations obtained from fits to electroweak data Note the correlation between 
the two masses in the latter. Within the current accuracy, the direct and indirect measure- 
ments are in approximate agreement. The central values of Mw and their errors, determined 
in several ways from indirect electroweak fits, are given in Table |I|. The results are evidently 
somewhat sensitive to the inclusion (or not) of data on the Z partial width ratios R;, and Rc 
and the SLD/SLC measurement of Am, all of which differ by 2.5 standard deviations or more 
from the Standard Model values. However, the conclusion on the agreement of the direct and 
indirect determinations is unchanged. As we shall see in the following sections, a significant 
reduction in the error on Mw is expected from both LEP2 and the Tevatron. 
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Figure 2: A comparison of direct versus indirect determinations of Mt and Mw- The contour 
for the indirect determination corresponds to a 70% confidence level. The dark shaded region 
within the contour is that compatible with the direct Higgs search limit, Mh > 65 GeV. 

1.3 Improved precision on Mw from the Tevatron 



The Tevatron data so far analysed, and shown in Fig. |I|, come from the 1992/3 data-taking 
(Run la). The resuhs from CDF [0] are based on approximately 19 pb~^ and are final, whereas 
those from DO [^] are based on approximately 13 pb~^ and are still preliminary. It is to be 
expected that the final result will have a smaller error. In addition, there will be a significantly 
larger data sample from the 1994/6 data-taking (Run lb). This should amount to more than 
100 pb~^ of useful data for each experiment. When these data are analysed it is envisaged that 
the total combined error on will be reduced to about 70 MeV. In particular, the combined 
CDF/DO result will depend on the common systematic error arising from uncertainties in the 
parton distribution functions. Thus when the first Mw measurements emerge from LEP2 one 
may assume that the world average error will have approximately this value. For more details 
see Ref. 
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After 1996 there will be a significant break in the Tevatron programme. Data-taking will 
start again in 1999 with a much higher luminosity (due to the main injector and other im- 
provements). Estimating the error on Mw which will ultimately be achievable (with several 
fb~^ of total luminosity) is clearly more difficult. If one assumes that an increase in the size of 
the data sample leads to a steady reduction in the systematic errors, one might optimistically 
envisage that the combined precision from the Tevatron experiments will eventually be in the 
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Table 1: Results of Standard Model fits to Mt, Mh and as{Mz)- The fits are to LEPl and SLC 
data. An upper limit Mh = 1000 GeV has been imposed. The electromagnetic coupling constant 
a(Mz) is also determined in these fits. Also given is the result of the fit given in terms of M^f^. 
The results are given for all data, as well as for excluding Rf, and Rc and also Alr. 



AM-w = 30 — 40 MeV range, assuming a common systematic error of about 25 MeV |T2 . 
However it is important to remember that these improved values will be obtained after the 
LEP2 measurements. 



1.4 Impact of a precision measurement of Mw 

Within the Standard Model, the value of Mw is sensitive to both Mt and Mh- For example, for 
a fixed value of Mh = 300 GeV, a precision of AM-w = ±25 MeV translates to a precision on 
Mt of ±3.9 GeV. The impact of a precise measurement of M-w on the indirect determination 
of Mh is shown in Table 





Mt = 180 GeV (fixed) 


Mt = 180 + 5 GeV 


Mh (GeV) 


AMw = 25 MeV 


AMw = 50 MeV 


AMw = 25 MeV 


AMw = 50 MeV 


100 


+48, -36 


+ 112, -63 


+86, -54 


+ 140, -72 


300 


+ 111, -84 


+259, -148 


+ 196, -126 


+323, -168 


800 


+297, -212 


+740, -367 


+538, -310 


+958, -413 



Table 2: Estimated error on Mh (in GeV) for several values 0/ Mh and for AMw = 25 and 
50 MeV. The estimates are given for both Mt = 180 GeV (fixed) and for Mt = 180 ± 5 GeV. 
All other Standard Model parameters are fixed. 



In order to assess the impact of a precise measurement of Mw it is necessary to make an 
estimate of the improvements which will be made on the electroweak data from LEPl and SLC. 
Details of the improvements which are assumed here are discussed in |]13[. The importance of 
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a precise measurement of Mw can perhaps best be appreciated by considering the (almost) 
model independent e parameters [|1^. The parameter ei (= Ap) is sensitive mainly to the Z 



partial and total widths. The parameter €3 depends linearly on both Ap and Ak, where Ak is 
determined from sin^ ^e//- The parameter €2 depends linearly on Ap, A/t and Ar^y. This latter 
quantity is determined essentially from Mw, and so improvements in the precision of €2 depend 
directly on improving the error on Mw- This is illustrated in Fig. ^, which shows the 70% 
confidence level contours for fits to projected global electroweak data. The different contours 
correspond to different values of AMw- In these fits all electroweak data measurements have 
been set to correspond to the Standard Model values Mt = 180 GeV, Mh = 100 GeV and 
as(Mz) = 0.123. The e variables are constructed to be sensitive to vector boson propagator 
effects, from both physics within the Standard Model and beyond. 




Figure 3: 70% confidence level contour plots for fits to t's in units of 10 ^ . The outer contour 
is for a fit with the current error, AMw = ±160 MeV, whereas the inner contours are for 
AMw = ±50 MeV and AMw = ±25 MeV respectively. 

Numerically, the projected data give a precision 

Aei = ±1.1 X 10~^ Aeg = ±1.0 x 10"^ . (3) 

For AMw = 160 MeV, the error Ae2 = ±3.7 x 10^'^ is obtained, whereas for the projected 
errors on Mw one obtains 

Ae2 = ±1.0 (±1.6) X 10~^ for AMw = 25 (50) MeV . (4) 

The smaller the volume in e space allowed by the precision electroweak measurements, the 
greater the constraint on physics beyond the Standard Model. 

The MSSM is arguably the most promising new-physics candidate. It is therefore especially 
important to consider the MSSM prediction for Mw. Figure ^ |l^ shows Mw as a function of 
Mt in the SM (solid lines) and in the MSSM (dashed lines). In each case the prediction is a 
band of values, corresponding to a variation of the model parameters (dominantly Mh in the 
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SM case, with 90 GeV < Mh < 1000 GeV chosen here) consistent with current measurements 
and hmits. An additional constraint of 'no SUSY particles at LEP2' is imposed in the MSSM 
calculation. 
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Figure 4: Predictions for My/ as a function of Mt in the SM (solid lines) and in the MSSM 
(dashed lines), from Ref. ]7^ /. In each case the prediction is a band of values, corresponding 
to a variation of the model parameters consistent with current measurements and limits. An 
additional constraint of 'no SUSY particles at LEP2' is imposed in the MSSM calculation. 



1.5 Methods for measuring 

Precise measurements of Mw can in principle be obtained using the enhanced statistical power 
of the rapidly varying total cross-section at threshold, the sharp (Breit-Wigner) peaking be- 
haviour of the invariant-mass distribution of the decay products and the sharp end-point 
spectrum of the lepton energy in decay. One can obtain a rough idea of the relative power 
of these methods by estimating their statistical precision assuming 100% efficiency, perfect 
detectors and no background. More complete discussions are given in Sections ^ and |^. 

A) Threshold cross-section measurement of the process e"'"e~ W"'"W~. The statis- 
tical power of this method, assuming 100% signal efficiency and no background, is 

AMw > 91 MeV^ ^°°^^ - , (5) 

where the minimum value is attained at y/s ~ 161 GeV. Here C denotes the total 
integrated luminosity. 
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B) Direct reconstruction methods, which reconstruct the Breit-Wigner resonant shape 
from the final states using kinematic fitting techniques to improve the mass resolution. 
The statistical power of this method, again assuming 100% efficiency, perfect detector 
resolution and no background, can be estimated as 

AMw~^-50Mev/i°^, (6) 

approximately independent of the collider energy. This order of magnitude estimate is 
confirmed by more detailed studies, see below. 

C) Determination of Mw from the lepton end-point energy. The end-points of the lepton 
spectrum in W^W~ qq/z/ depend quite sensitively on the W mass. For on-shell W 
bosons at leading order: 

E^<Ei<E+, E± = ^ (^1± ^1 - 4M2,/s) . (7) 

In this case the statistical error on Mw is determined by the statistical error on the 
measurement of the lepton end-point energy, 

Js - 4M2 

AMw = ^ AE± . (8) 

In practice, however, the end-points of the distribution are considerably smeared by finite 
width effects and by initial state radiation, and only a fraction of events close to the 
end-points are sensitive to Mw- This significantly weakens the statistical power of this 
method from what the naive estimate (|) would predict. 

The detailed studies described in the following sections show that the errors which can real- 
istically be achieved in practice are somewhat larger than the above estimates for Methods 
A and B. The statistical precisions of the two methods are in fact more comparable (for the 
same integrated luminosity) than the factor 2 difference suggested by the naive estimates (^) 
and (^). The overall statistical error for Method C has been estimated at (9(300 MeV) |1[] for 
C = 500 pb~^ , significantly larger than that of the other two methods. It will not therefore 
be considered further here, although it is still a valid measurement for cross-checking the other 
results. 

It is envisaged that most of the LEP2 data will be collected at energies well above threshold, 
and so the statistically most precise determination of Mw will come from Method B. However 
with a relatively modest amount of luminosity spent at the W"*"W~ threshold (for example 
50 pb~^ per experiment). Method A can provide a statistical error of order 100 MeV, not 
significantly worse than Method B and with very different systematics. The two methods can 
therefore be regarded as complementary tools, and both should be used to provide an internal 
cross-check on the measurement of the W mass at LEP2. This constitutes the main motivation 
for spending some luminosity in the threshold region. 
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The threshold cross-section method is also of interest because it appears to fit very well into 
the expected schedule for LEP2 operation in 1996. It is anticipated that the maximum beam 
energy at LEP2 will increase in steps, with the progressive installation of more superconducting 
RF cavities, in such a way that a centre-of-mass energy of 161 GeV will indeed be achievable 
during the first running period of 1996. This would then be the ideal time to perform such a 
threshold measurement. The achievable statistical error on Mw depends of course critically on 
the available luminosity at the threshold energy. In Section ^ we present quantitative estimates 
based on integrated luminosities of 25, 50 and 100 pb~^ per experiment. 

1.6 Theoretical input information 

1.6.1 Cross-sections for the W+W signal and backgrounds 

Methods (A) and (B) for measuring Mw described above require rather different theoretical 
input. The threshold method relies on the comparison of an absolute cross-section measurement 
with a theoretical calculation which has Mw as a free parameter. The smallness of the cross- 
section near threshold is compensated by the enhanced sensitivity to Mw in this region. In 
contrast, the direct reconstruction method makes use of the large W"'"W~ statistics at the higher 
LEP2 energies, ^175 GeV. Here the more important issue is the accurate modeling of the 
line-shape, i.e. the distribution in the invariant mass of the decay products. 

In this section we describe some of the important features of the theoretical cross-sections 
which are relevant for the M^ measurements. A more complete discussion can be found in the 
contribution of the WW and Event Generators Working Group to this Report [|T^] . 

We begin by writing the cross-section for e"'"e~ — 4f (+7,g, . . .), schematically, as 

aww = 0-^^ (1 + ^Ew + 5qcd) , (9) 

We note that this decomposition of the cross-section into 'signal' and 'background' contribu- 
tions is practical rather than theoretically rigorous, since neither contribution is separately 
exactly gauge invariant nor experimentally distinguishable in general. The various terms in (^) 
correspond to 

(i) 0"^^: the Born contribution from the 3 'CC03' leading-order diagrams for e+e^ — > 
W'^W" involving t-channel u exchange and s-channel 7 and exchange, calculated 
using off-shell W propagators. 

(ii) 5ew- higher-order electroweak radiative corrections, including loop corrections, real pho- 
ton emission, etc. 

(iii) Sqcd- higher-order QCD corrections to W+W~ final states containing qq pairs. For 
the threshold measurement, where in principle only the total cross-section is of primary 
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interest, the effect of these is to generate small corrections to the hadronic branching ratios 
which are entirely straightforward to calculate and take into account. More generally, such 
QCD corrections can lead to additional jets in the final state, e.g. W"'"W~ qqqqg from 
one hard gluon emission. This affects the direct reconstruction method, insofar as the 
kinematic fits to reconstruct M-w assume a four-jet final state, and both methods insofar 
as cuts have to be imposed in order to suppress the QCD background (see Sections 
below). Perturbative QCD corrections, real gluon emission to O(a^) and real plus virtual 
emission to 0{as), have been recently discussed in Refs. |TP|, ^ respectively, together 
with their impact on the measurement of Mw. 

(iv) (Tbkd: 'background' contributions, for example from non-resonant diagrams (e.g. e'^e~ — >■ 
/i~^t'^W~) and QCD contributions e+e" qqgg(7), qqqq(7) to the four-jet final state. 
All of the important backgrounds have been calculated, see Table ^ below. At threshold, 
the QCD four-jet background is particularly large in comparison to the signal. 

In what follows we consider (i) and (ii) in some detail. Background contributions and how to 
suppress them are considered in later sections. 



1.6.2 The W+W ofF-shell cross-section 



The leading-order cross-section for off-shell W W production was first presented in Ref. |2l 

a{s) = dSi dS2 p{si)p{s2) (To{s,Si,S2) , (10) 

JO Jo 

where 

( \ = }lIjL '1 (-]-]] 

The cross-section ctq can be written in terms of the z/, 7 and Z exchange contributions and 
their interferences: 

cro(s. Si, S2) = 7—7 — 5 [a^^ + azz + a-iZ + a-w + ^uz + 01.7] , (12) 

zbbvrs S\S2 

where = e^/sin^ 9w- Explicit expressions for the various contributions can be found in 
Ref. |21| for example. The stable (on-shell) W^W^ cross-section is simply 



a-(s) = oro(s,M^,M^). (13) 

A theoretical ansatz of this kind will be the basis of any experimental determination of the 
mass and width of the W boson. The reason for this is the large effect of the virtuality of the 
W bosons produced around the nominal threshold. An immediate conclusion from Eq. ([I0|) 
may be drawn: the W mass influences the cross sections exclusively through the off-shell W 
propagators; all the other parts are independent of Mw and Fw (neglecting for the moment 
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the relatively minor dependence due to radiative corrections) . It will be an important factor in 
the discussion which follows that near threshold the (unpolarized) cross-section is completely 
dominated by the t-channel neutrino exchange diagram. This leads to an S-wave threshold 
behaviour at ~ /3, whereas the s-channel vector boson exchange diagrams give the characteristic 
P-wave behaviour cr^ ~ 

By tradition (for example at LEPl with the Z boson), the virtual W propagator in Eq. ([Tl|) 
uses an s-dependent width, 

rw(5) = Tw , (14) 

where Fw = rw(^w)- Another choice, equally well justified from a theoretical point of view, 
would be to use a constant width in the W propagator (for a discussion see Ref. [P^): 



W-L w 



(15) 



The numerical values of the width and mass in the two expressions are related p^: 



1 f2 

Mw = Mw - = Mw - 26.9 MeV , (16) 

Tw = rw-^]J = rw-0.7MeV. (17) 

These relations may be derived from the following identity: (s — m\j + iMwEw) = (s — 
+ isVy^ / My^) / {1 + iVyj/M-^). Numerically, the consequences are below the anticipated 
experimental accuracy. 



1.6.3 Higher-order electroweak corrections 

The complete set of 0{a) next-to- leading order corrections to W^W~ production has been 



calculated by several groups ||2J, for the on-shell case only, see for example Refs. ||26| , |TS 
and references therein. There has been some progress with the off-shell (i.e. four fermion 
production) corrections but the calculation is not yet complete. However using the on-shell 
calculations as a guide, it is already possible to predict some of the largest effects. For ex- 
ample, it has been shown that close to threshold the dominant contribution comes from the 
Coulomb correction, i.e. the long-range electromagnetic interaction between almost stationary 
heavy particles. Also important is the emission of photons coUinear with the initial state e^ 
('initial state radiation') which gives rise to logarithmic corrections ~ a\n{s/rn^). These lead- 
ing logarithms can be resummed to all orders, and incorporated for example using a 'structure 
function' formalism. In this case, the generalization from on-shell to off-shell W^s appears to be 
straightforward. For the Coulomb corrections, however, one has to be much more careful, since 
in this case the inclusion of the finite W decay width has a dramatic effect. Finally, one can 
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incorporate certain important higher-order fermion and boson loop corrections by a judicious 
choice of electroweak couphng constant. Each of these effects will be discussed in turn below. 

In summary, certain 0{a) corrections are already known to be large because their coefficients 
involve large factors like \n{s/ml), ^Mw/Fw, /My^, etc. Once these are taken into account, 
one can expect that the remaining corrections are no larger than 0{a). When estimating the 
theoretical systematic uncertainty on the W mass in Section ^ below, we shall therefore assume 
a conservative overall uncertainty on the cross-section of ±2% from the as yet uncalculated 
0{a) and higher-order corrections. 



1.6.4 Coulomb corrections 

The Coulomb corrections for on-shell and off-shell W"'"W~ production have been discussed in 
detail in Refs. [^, |2^, |2^, where a complete set of references to earlier studies can also be 
found. 



The result for on-shell W^W~ production is well-known [^] — the 0{a) correction diverges 
as aw/vo as the relative velocity Vq = 2[1 — 4M-^/s]^/^ of the W bosons approaches zero at 
threshold. Note that cxo ~ Vq near threshold and so the Coulomb-corrected cross-section is 
formally non- vanishing when y/s = 2Mw. 

For unstable W'^W" production the finite decay width F-w screens the Coulomb singularity 
^1, so that very close to threshold the perturbative expansion in air/vQ is effectively replaced 



by an expansion in aix^ M^^ /T\^ [^. In the calculations which follow we use the expressions for 
the 0{a) correction given in Ref. [^]. The net effect is a correction which reaches a maximum 
of approximately +6% in the threshold region. Although this does not appear to be large, we 
will see below that it changes the threshold cross-section by an amount equivalent to a shift in 



Mw of order 100 MeV. In Ref. the 0(a) result is generalized to all orders. However the 
contributions from second order and above change the cross-section by ^ 1% in the threshold 
region |^ (see also |^) and can therefore be safely neglected. 

Note also that the Coulomb correction to the off-shell W+W" cross-section 

cr(s) = / dsi ds2 p{si)p{s2) cro{s,Si,S2)[l + Sc{s,Si,S2)] (18) 

Jo Jo 

provides an example of a (QED) interconnection effect between the two W bosons: the exchange 
of a soft photon distorts the line shape {da /d^/si) of the and therefore, at least in principle. 



affects the direct reconstruction method [^, |33|, |3^. In Ref. |3^, for example, it is shown 
that the Coulomb interaction between the W bosons causes a downwards shift in the average 
reconstructed mass of 0{■naV^^) ~ 0(20 MeV). Selecting events close to the Breit-Wigner 
peak reduces the effect somewhat. However the calculations are not yet complete, in that QED 
interactions between the decay products of the two W bosons are not yet fully included. 
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1.6.5 Initial state radiation 



Another important class of electroweak radiative corrections comes from the emission of photons 
from the incoming e"*" and e~. In particular, the emission of virtual and soft real photons 
with energy E < uj gives rise to doubly logarithmic contributions ~ a\\i{s/m1)\n{s/uj'^) at 
each order in perturbation theory. The infra-red (Incu) logarithms cancel when hard photon 
contributions are added, and the remaining coUinear (ln(s/mg)) logarithms can be resummed 
and incorporated in the cross-section using a 'flux function' or a 'structure function' [^] (see 



also Refs. pq, 38, 39, 18 



The ISR corrected cross section in the flux function (FF) approach is 



dsids2 
where x = 1 — s' / s and 

with 



ds' { V , 

— I s)[cTcCn(s',Si,S2) + 5c CTccs (s'. Si , S2) 
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non— umv 



ISR 



Fix, s) 



t 



tx''\l + S) + H{s',s) 



2a 

IT 



In 



1 



(19) 



(20) 



(21) 



The S term comes from soft and virtual photon emission, the Ti term comes from hard photon 
emission, 6c contains the Coulomb correction (0), Uccs is the doubly resonating Born cross 
section, and cxccn the background contributions. Explicit expressions can be found in the above 
references. The additional term crfsjlj^""^ is discussed in Refs. ||38| and [|l8l . The invariant mass 
lost to photon radiation may be calculated as 



K) = IJ dsids, J (l - - 



da 



s J dsids2ds' ' 

where da / dsids2ds' is the contents of the curly brackets in Eq. (|19|). 
Alternatively, the structure function (SF) approach may be used: 

daQ-Enis) 



[22] 



dsids2 



J dxi J dX2D{xi,s)D{x2,s)<acCniXiX2S,Si,S2)+ Sc CTCCS} 

-min„min ^ ' 



where 



D(x, s) 



-(l-xy/'-\l + S) + Hix,s). 



Here, the invariant mass loss is 



^ f ^ A f A A nf \T^f ^^/^n da(xiX2S, 81,82) 

- / d8id82 / dXidX2D{Xi,8)D[X2,8) — [l - X1X2) — ; ; ; ; 

a J J 2 d8id82dxidx2 



(23) 



(24) 



(25) 
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In addition, the radiative energy loss may be determined, 



1 /" f 

{E^) = —J dsids2 J dxidx2 D{xi, s)D{x2, s)-^{2 — xi 



da 



- X2, 



dsids2dxidx2 



(26) 



Initial state radiation affects the W mass measurement in two ways. Close to threshold the 
cross-section is smeared out, thus reducing the sensitivity to Mw (see Fig. ^ below). For the 
direct reconstruction method, the relatively large average energy carried away by the radiated 
photons leads to a large positive mass-shift if it is not taken into account in the rescaling of 
the final-state momenta to the beam energy (see Section ^ below) . By rescaling to the nominal 
beam energy we obtain for the mass-shift AMw = {E^)M^^/^/s. Note however that a fit to the 
mass distribution gives more weight to the peak, and therefore in practice the effective value 
of {E^) or (m^) is less than that given by Eqs. (^,|^,^) (see Section [3l^ ). Table ^ shows 
the influence of the various cross-section contributions on the average energy and invariant 
mass losses. The invariant mass loss may be calculated both in the SF and FF approaches. A 
comparison shows that the predictions in both schemes differ only slightly, which allows us to 
use the numerically faster FF approach for the numerical estimates. At the lower LEP2 collider 
energies, the energy and invariant mass losses are nearly equal, while at higher energies their 
difference is non-negligible. Note also that the inclusion of the non-universal ISR corrections 
and background terms is of minor influence. The latter has been studied only for CCll processes; 
for reactions of the CC20 type the background is larger and the numerical estimates are not 
yet available. The Coulomb correction is numerically important and cannot be neglected |29[ . 



The dependence of the predictions on the details of the treatment of QED is discussed in detail 
in Ref. |jl8| and will not be repeated here. 





(GeV) 


175 


192 


205 


SF, CCS 


a (pb) 


13.182 


16.488 


17.077 




{E,) 


1115 


2280 


3202 






1112 


2271 


3185 






3 


9 


17 


change to FF 


A(m^) 


0.5 


-0.8 


-4.2 


add 6c 


A(m^) 


11.8 


16.3 


20.8 


add disR 


A(m^) 


0.2 


0.3 


0.3 


add (Tccn 


A(m^) 


<0.1 


0.2 


0.5 



Table 3: Influence of different cross-section treatments on the average energy loss {E^) and 
invariant mass loss (m^) in MeV. 
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parameter 


value 


Mz 


91.1888 


Mw 


80.23 


Tz 


2.4974 


Fw 


2.078 




137.0359895 




1.16639 X 10-5 GeV-2 


Sin y^y = sm (7^ 


0.2320 


me 


5.1099906 X 10-4 


(^c)2 


3.8937966 x 10^ pb GeV^ 



Table 4: Parameter values used in the numerical calculations in Table and Fig. ^. Masses 
and widths are given in GeV. 

1.6.6 Improved Born approximation 



In the Standard Model, three parameters are sufficient to parametrise the electroweak interac- 
tions, and the conventional choice is {a, G^, Mz} since these are the three which are measured 
most accurately. In this case the value of Mw is a prediction of the model. Radiative corrections 
to the expression for Mw in terms of these parameters introduce non-trivial dependences on 
Mi and Mh, and so a measurement of Mw provides a constraint on these masses. However 
the choice {a,G^,Mz} does not appear to be well suited to W"'"W- production. The reason 
is that a variation of the parameter Mw, which appears explicitly in the phase space and in 
the matrix element, has to be accompanied by an adjustment of the charged and neutral weak 
couplings. Beyond leading order this is a complicated procedure. 



It has been argued |^0| that a more appropriate choice of parameters for LEP2 is the set 
{Mw,G^,Mz} (the so-called G^-scheme), since in this case the quantity of prime interest is 
one of the parameters of the model. Using the tree-level relation 



9 



V sin^ d 



w 



aV2G^MI 



w 



(27) 



we see that the dominant t-channel neutrino exchange amplitude, and hence the corresponding 
contribution to the cross-section, depends only on the parameters Mw and G^. It has also been 
shown 1^ that in the G^-scheme there are no large next-to-leading order contributions to the 
cross-section which depend on Mt, either quadratically or logarithmically. One can go further 
and choose the couplings which appear in the other terms in the Born cross-section such that 
all large corrections at next-to-leading order are absorbed, see for example Ref. |4l[]. However 
for the threshold cross-section, which is dominated by the t-channel exchange amplitude, one 
can simply use combinations of and defined by Eq. (E^) for the neutral and charged weak 
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couplings which appear in the Born cross-section, Eq. (|T2[). 

In summary, the most model-independent approach when defining the parameters for com- 
puting the e+e~ — > W"'"W~ cross-section appears to be the G^-scheme, in which Mw appears 
explicitly as a parameter of the model. Although this makes a non-negligible difference when 
calculating the Born cross-section, compared to using a and sin^ 6^/ to define the weak cou- 
plings (see Table § below), a full next-to-leading-order calculation will remove much of this 
scheme dependence 





v/i = 161 GeV 


v/i = 175 GeV 


(Jo (on- shell, a) 


3.813 


15.092 


(To (on- shell, G^) 


4.402 


17.425 


ao (off-shell with rw(M^), G^) 


4.747 


15.873 


(Jo (off-shell with rw(s), G^) 


4.823 


15.882 


. . . + 0{a) Coulomb 


5.122 


16.311 


... + 0(a) Coulomb + ISR 


3.666 


13.620 



Table 5: Decomposition of the theoretical e~^e —>■ W^W cross-section (in picobarns) as 
defined and discussed in the text, at two LEP2 collider energies. 



1.6.7 Numerical evaluation of the cross-section 

Figure |^ shows the e"'"e~ — > W"^W~ cross-section at LEP2 energies. The different curves 
correspond to the sequential inclusion of the different effects discussed above. The parameters 
used in the calculation are listed in Table ^. Note that both the initial state radiation and the 
finite width smear the sharp threshold behaviour at = 2Mw of the on-shell cross-section. 
The different contributions are quantified in Table |^, which gives the values of the cross-section 
in different approximations just above threshold {^/s = 161 GeV) and at the standard LEP2 
energy of -y/i = 175 GeV. At threshold we see that the effects of initial state radiation and 
the finite W width are large and comparable in magnitude. For the threshold method, the 
primary interest is the dependence of the cross-section on Mw. This will be quantified in 
Section ^ below. For both methods, the size of the background cross-sections is important. For 
completeness, therefore, we list in Table ^ some relevant cross-section values obtained using the 
PYTHIA Monte Carlo. This includes finite-width effects, initial state radiation and Coulomb 
corrections. Notice that the values for aww agree to within about 1% accuracy with those 
given in the last row of Table ^. 
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Reaction 


Cross-section (pb) 




at 161 GeV 


at 175 GeV 


at 192 GeV 


e+e- ^ W+W- ^ all 


3.64 


13.77 


17.10 


e+e" — > W+W"— >qqqq 


1.67 


6.30 


7.85 


e+e" — >• W+W-^qq/z/ 


1.59 


6.02 


7.46 


c+c- ^ W+W-^/z//z/ 


0.38 


1.44 


1.79 


e+e^ — >• Z/7 — >• all 


221. 


172. 


135. 


e+e" — > Z/7 — > qq 


151. 


116. 


91. 


e+e- ^ Z/7^ Z+r,i/I7 


70. 


56. 


44. 


e+e- ^ (Z/7) (Z/7) ^ all 


0.46 


0.44 


1.12 


e+e- Zee ^ all 


2.53 


2.70 


2.85 


e+e" Wcu all 


0.37 


0.51 


0.72 


e+e^ — s> e+e^, cos 6' < 0.95 


351. 


297. 


247. 



Table 6: Cross-sections (in picobarns) for signal and background channels as given by PYTHIA. 
W+W", (Z/7) (Z/7), Wez/ and Zee refer to four-fermion production with intermediate formation 
of either one or two vector bosons; in the case of Z pair production, the generator includes the 
photon contribution. The process Z/7 all refers to the production of a fermion pair via a 
Z boson or photon. All interactions include initial state radiation, and therefore may contain 
additional photons in the final state. The final row lists cross-sections generated for the Bhabha 
scattering process (including t-channel exchange) using the BABAMC program. 
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Figure 5: The cross-section for e^e'^ — > W^W~ in various approximations: (i) Born (on-shell) 
cross-section, (ii) Born (off-shell) cross-section, (Hi) with first order Coulomb corrections, and 
(iv) with initial state radiation. The parameter values are listed in Table 0. 

2 Measurement of from the W+W~ Threshold Cross- 
Sectionf] 



As discussed in Section [L^, one can exploit the rapid increase of the W^'^W^ production cross- 
section at y/s ~ 2Mw to measure the W mass. In the following, we briefly discuss the basic 
features of this method, suggest an optimal collider strategy for data-taking, and estimate the 
statistical and systematic errors. The intrinsic statistical limit to the resolution on Mw is 
shown to be energy- dependent: in particular, arguments are presented in favour of a single 
cross-section measurement at a fixed energy -^i ~ 161 GeV. 



^prepared by D. Gele, T.G. Shears, W.J. Stirling, A. Valassi, M.F. Watson 
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2.1 Collider strategy 



The cross-section for W'^'W^ production increases very rapidly near the nominal kinematic 
threshold y/s = 2Mw, although the finite W width and ISR smear out the abrupt rise of the 
Born on-shell cross-section. This means that for a given a/s near threshold, the value of the 
cross-section is very sensitive to Mw- This is illustrated in Fig. where the W^W~ excitation 
curve is plotted for various values of the W mass. The calculation is the same as that discussed 
in Section |L6| , and includes finite W width effects, ISR and QED Coulomb corrections. A 
measurement of the cross-section in this region therefore directly yields a measurement of Mw- 




Figure 6: The cross-section for W~^W~ production as a function of ^/s in the threshold region, 
for various values of Myj. Finite-width effects, the QED Coulomb correction and Initial State 
Radiation are included. 

For an integrated luminosity C and an overall signal efficiency eww = CjERj (where the 
sum extends over the various channels selected, with branching ratios BRj and efficiencies e^), 
the error on the W"'"W~ cross-section due to signal statistics is given by 



Aa„,„. = = (28) 
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where N = eww^ww^ is the number of selected signal events. The corresponding error on the 
W mass is 

Ailf-W = y/O'wW 

The sensitivity factor ^yaww\dM^fil / daww\ is plotted in Fig. ^ as a function of ^/s — 2M^^. 
There is a minimum at 

(v/i) ~ 2Mw + 0.5 GeV, (30) 

corresponding to a minimum value of approximately 0.91 GeV pb~^^^. Note that the 0.5 GeV 
offset of the minimum of the sensitivity above the nominal threshold is insensitive to the actual 
value of Mw, since in the threshold region the cross- section is to a first approximation a function 
of a/s — 2Mw only. 

As discussed below, the statistical uncertainty is expected to be the most important source 
of error for the threshold measurement of Mw: the optimal strategy for data-taking consists 
therefore in operating at the collider energy (a/s) in order to minimize the statistical error 
on Mw. The statistical sensitivity factor is essentially flat within (y/s) ± 0.65 GeV, where it 
mcreases at most to 0.95 GeV ph'^^^ (+4%); bearing in mind that the present uncertainty on 
Mw from direct measurements is 160 MeV (and is expected to decrease further in the coming 
years), this corresponds to ±2a on the current world average Mw value. In other words, Mw is 
already known to a level of precision good enough to choose, a priori, one optimal energy for the 
measurement of the W"'"W~ cross-section at the threshold. Using the latest world average value 
Mw = 80.26 ±0.16 GeV (see Eq. (|lD) gives an optimal colhder energy of {^/s)°P^ ~ 161.0 GeV. 



w 



da 



WW 



V^ww^ 



(29) 



2.2 Event selection and statistical errors 

The error on the W mass due to the statistics of W^W~ events collected has been given in the 
previous section. Background contamination with an effective cross-section cTbkg introduces an 
additional statistical error. The overall effect is that the statistical error on Mw is modified 
according to 

AMw ^ AMw Jl + ""^^^ . (31) 

V ^WW^WW 

In the following subsections we present estimates of this statistical error for realistic event 
selections, for an integrated luminosity of 100 pb~^ at ^/s = 161 GeV. Tight selection cuts 
are required to reduce the background contamination while retaining a high efficiency for the 
signal, especially as the signal cross-section is a factor of 4-5 lower at threshold than at higher 
centre-of-mass energies. 

The studies are based on samples of signal and background events generated by means 
of Monte Carlo programs (mainly PYTHIA [^) tuned to LEPl data. These events were 
run through the complete simulation program giving a realistic detector response, and passed 
through the full reconstruction code for the pattern recognition. 



22 




Figure 7: The sensitivity of the W^W~ cross-section to the W mass, plotted as a function of 
^/s — 2Mw. The significance of the three curves to the W mass measurement is discussed in 
the text. A value of M^^ = 80.26 GeV has been used in the calculations. 

2.2.1 Fully hadronic channel, W"'"W~— >qqqq. 

The pure four quark W"'"W~ decay mode benefits from a substantial branching ratio (46%) 
corresponding to a cross-section ct(W+W~— *qqqq; y/s = 161 GeV) = 1.67 pb. Obviously, the 
typical topology of such events consists of four or more energetic jets in the final state. Due 
to its large cross-section (see Table |), the main natural background to this four-jet topology 
comes from e"'"e~ — >• Z(— qq(+ng))7 events which can be separated into two classes depending 
on the virtuality of the Z: (i) the production of an on- mass-shell Z accompanied by a radiative 
photon of nearly 55 GeV (at s/s = 161 GeV), which is experimentally characterised by missing 
momentum carried by the photon escaping inside the beam pipe (typically 70% of the time), 
and (ii) events with a soft ISR 7 and a large total visible energy, which potentially constitute 
the most dangerous QCD background contribution. 

Note that a semi-analytical calculation of the genuine four-fermion background cross-section 
[|T6| for a wide range of four-fermion final states (with non-identical fermions) shows that in 



the threshold region c'"bkgd(4f)/crvi/vv' 1%; and therefore the effect on the Mw determination 
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from these final states is negligible. 

Although the effective four-jet-like event selection depends somewhat on the specific detector 
under consideration, a general and realistic guideline selection can be described. The most 
relevant conditions to be fulfilled by the selected events can be summarised as follows: 

• A minimum number of reconstructed tracks of charged and neutral particles is required. 
A typical value is 15. This cut removes nearly all low multiplicity reactions, such as 
dilepton production (e"'"e~ —>•//,/ = e, /x, r) and two-photon processes. 

• A veto criterion against hard ISR photons from qq7 events in the detector acceptance 
can be implemented by rejecting events with an isolated cluster with significant electro- 
magnetic energy (larger than 10 GeV for example). 

• A large visible energy, estimated using the information from tracks of charged particles 
and from the electromagnetic and hadron calorimeters. For example, a minimum energy 
cut value of 130 GeV reduces by a factor of 2 the number of qq7 events with a photon 
collinear to the beam axis. 

• A minimum number (typically 5) of reconstructed tracks per jet. This criterion acts on 
the low multiplicity jets from r decays as well as from 7 conversions or 7 interactions 
with the detector material. 

• A minimum jet polar angle. The actual cut value depends on the detector setup, but is 
likely to be around 10° — 15°. This cut is mainly needed to eliminate poorly measured jets 
in the very forward region, where the experiments are generally less well instrumented. 

These selection criteria almost completely remove the harmless background sources (e"'"e~ —>■ 
ZZ, Zee and Wei/), but there is still an unacceptable level of e+e" Z(^ <iq(+^g))7 contami- 
nation (about three times higher than the signal). The second step is to suppress the remaining 
QCD background (from e+e~ — > o^gg'^soft events) by performing a mass reconstruction 
based on a constrained kinematic fit. The following additional criteria can then be imposed: 

• A probability cut associated with a minimum constrained dijet mass requirement - 
a typical choice of standard values of 1% and 70 GeV respectively is used here. This 
procedure appears to be an efficient tool to improve the mass resolution and therefore to 
reduce the final background, see Fig. |^. 

In summary, a reasonable signal detection efficiency in excess of 50% is achievable for 
W"''W~— >qqqq events. Although the final rejection factor of QCD events is approximately 
500, a substantial residual four-jet background still remains, giving a purity of around 70%. 
The contributions from other backgrounds (ZZ, Zee, Wez/ and two-photon events) are negligible 
for the four-jet analysis. 
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Figure 8: Distributions of the fitted dijet mass using a constrained kinematic fit as described 
in the text. The solid histogram shows the W^W~ — > qqqq signal and the hatched histogram 
represents the e"^e~ —>■ qq'j background events. 



The signal and background efficiencies for the typical event selection described above are 
given in Table ^ assuming a total integrated luminosity of 100 pb~^ (i.e. 25 pb~^ per interaction 
point). 





W+W ^qqqq 


W+W-^qq/z/ 




Signal cross-section 
Signal efficiency 
AMw (stat.) for signal 


0.94 pb 

55% 
180 MeV 


0.76 pb 

47% 
197 MeV 


0.23 pb 

60% 
354 MeV 


Background cross-section 
AMw (stat.) for background 


0.39 pb 
106 MeV 


0.03 pb 
37 MeV 


0.01 pb 

74 MeV 


Total AMw (stat.) 


209 MeV 


200 MeV 


360 MeV 



Table 7: Typical accepted cross-sections at y/s = 161 GeV and corresponding statistical uncer- 
tainty on Mw, assuming an integrated luminosity of 100 pb~^ . Note that I = e, /x, r. 



2.2.2 Semileptonic channel, W"'"W ^qqlu. 

The decay channel W"'"W~— >qq/i/ is characterised by the presence of two or more hadronic 
jets, an isolated, energetic lepton or a narrow jet in the case of hadronic r decays, and missing 
energy and momentum due to the undetected neutrino. Since the W"'"W~ cross-section is small 



25 



at threshold, processes such as e"'"e~ -^"Lj^ and e"'"e~ Zee constitute important backgrounds. 
Their cross-sections at 161 GeV are hsted in Table |^. Typical experimental selection cuts for 
the W"''W~^qqez/ and W"''W~^qq/iz/ channels include: 

• A multiplicity cut to reject purely leptonic events. By requiring at least 6 charged tracks 
in the event, backgrounds from the W"^W~^/z//z/, Z/7 — > and leptonic two-photon 
channels can be removed. 

• Identification of an electron or muon using standard experimental cuts. The lepton is 
required to have a high momentum and to be isolated from the hadronic jets. This 
isolation can be achieved by restricting the energy in a cone around the lepton candidate, 
or by requiring a minimum transverse momentum relative to the nearest jet. An example 
of such a cut is to require less than 2.5 GeV of energy (charged and neutral) inside a 
200 mrad cone around the track. This suppresses hadronic background, much of which 
originates from the semi-leptonic decays of heavy quarks inside jets. 

• Cuts on missing momentum and its direction. The neutrino carries away momentum, 
leaving the event with a net momentum imbalance. In order to distinguish the signal 
events from TP /'^ or two photon backgrounds, in which the missing momentum tends 
to be along the beam direction, the event is required either to have a large transverse 
momentum, or to have missing momentum pointing away from the beam direction. These 
quantities are illustrated in Fig. |^(a) and (b) and show clear differences between signal and 
background events. Effective experimental cuts are | cos 6*11118 1 < 0.95 and Pmis/-E'cm > 0.07. 

• Cuts on the corrected, visible energy in the event. Requiring i?vis/-E'cm < 1 removes some 
residual background from Z/7 — >■ hadrons. Similar results can be achieved by cutting on 
the energy or mass of the hadronic system. 

Additional improvements to the selection may be obtained by applying a kinematic fit, in 
which energy-momentum constraints are applied in conjunction with assumptions about the W 
masses. The background events tend to give poor probabilities in the fit. 

Selection cuts based on these quantities give efficiencies of 70-75% for the qqei/ and qq/iz/ 
channels, but only about 5% for the qqrz/ decays. The purity of the selected sample is 90- 
95%, with typical accepted cross-sections as listed in Table ^ The corresponding statistical 
uncertainty on Mw is approximately 200 MeV for 100 pb~^. More work will be needed to 
improve the efficiency for selecting r decays, and hence to enhance the statistical precision of 
the cross-section measurement. Hadronic r decays can be identified as a third jet with low 
multiplicity and a small opening angle. 

2.2.3 Fully leptonic channel, W+W--*/z//z/. 

The fully leptonic channel is not used by the direct reconstruction method due to the lack of 
sufficient constraints on the kinematics of the event, but it can be used for event-counting in 
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Icos 9.1 Transverse momentum (GeV) 



Figure 9: Examples of the momentum imbalance in candidate events for the W~^W~ ^qqHu 
channel: (a) direction of the missing momentum, (b) total momentum transverse to the beam 
direction. The open histogram indicates events which contain an isolated, identified lepton and 
pass all other selection cuts; the background from the Z/7, ZZ, Zee andWeu channels is shaded. 
All samples include detector simulation. 

the threshold measurement of Mw- It has, however, a branching ratio of only 11% (/ = e, /i, r), 
a factor 4 lower than the other two channels. This is reflected in the relative weight of this 
channel to the overall precision on Mw using the threshold method. 

The W^W^ ^lulu channel is characterised by two acoplanar, energetic leptons and large 
missing momentum. In 4/9 (1/9) of the cases, however, one (both) of these leptons is a r, which 
typically decays to a narrow hadronic jet. Typical experimental selections for all W^W^^lulu 
channels require: 

• Low charged multiplicity (typically not greater than 6), which allows the rejection of most 
of the hadronic backgrounds. The most important background to this channel is then 
given by dilepton (Z/7, two-photon or Bhabha) events. At least two good charged tracks 
(with a typical minimum energy of 1 GeV) are required in any case. 

• One identified electron or muon (requiring two would lose most of the events with one 
W — > TUr decay); this lepton is required to have an angle of typically at least 25° to 
the beam both Bhabhas and two-photon events are concentrated at low angles. 
Since the leptons from W+W~ decays at threshold have approximately half the beam 
energy, an energy window may be imposed on the lepton, centered around this value 
(such as [28-55] GeV); this is effective both against low-energy two-photon events (which 
can be further reduced by requiring a large total visible mass), and against di-muons and 
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Bhabhas, where the energy of each lepton is typically equal to the beam energy. 

• Explicit tagging of events with isolated, energetic photons or luminosity clusters allows 
one to reject radiative Z events with hard, detected ISR. 

• Cuts on the missing momentum and its direction can also be used: large transverse missing 
momentum (typically more than 4 GeV) is required to discriminate against two-photon 
events, while combined cuts on the angle between the jets and on the missing momentum 
out of the beam-thrust plane are very effective against di-tau Z/7 events. 

• To reconstruct the two original lepton directions, the event is forced into two jets, which 
must be of low invariant mass and well contained in the detector. A cut on the acopla- 
narity A0 between these two jets (such as A0 < 174°) is effective against all residual 
backgrounds. The distribution of this variable for signal and background events is shown 
in Fig. 0. 

Selection cuts of this kind give overall efficiencies of approximately 60% (60% to 70% in all 
individual channels except the tVt-tVj. channel, where it is only of the order of 30%), with a 
purity close to 95%. Typical accepted cross-sections are listed in Table 0. The corresponding 
statistical error on Mw is approximately 360 MeV (354 MeV from signal statistics plus 74 MeV 
from background statistics) for a total integrated luminosity of 100 pb^^. 




20 40 60 80 100 120 140 160 180 



A^(degrees) 

Figure 10: Distribution of the jet acoplanarity A0, in the selection of W~^W~ -^lulu events at 
threshold. The solid histogram includes the events selected by all cuts (except that on A(j)). The 
shaded area corresponds to the background from the Z/7 and 77 dilepton events. 
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2.3 Systematic errors 



Uncertainties in the W+W^ production cross-section translate directly into systematic errors 
on the W mass. The uncertainties fall into three categories: multiplicative uncertainties in 
the cross-section; an additive uncertainty due to background subtraction; other sources, for 
example beam energy and W width. 



2.3.1 Luminosity, higher-order corrections and selection efficiency 

Quantities which enter multiplicatively in the calculation or measurement of the cross-section 
contribute an error to the W mass which can be expressed as 



AM. 



w 



da 



dM. 



w 



a 



AC 



(32) 



where AC is the uncertainty on the multiplicative quantity C, and the cross-section a includes 
contributions from signal and background weighted by their relative efficiencies. The quantity 
a\da / dMy^l'^ is shown in Fig. as a function of ^/s. Like the statistical sensitivity factor, it 
also has a minimum near the nominal threshold and has the value of approximately 1.7 GeV 
at (ys)°P* ~ 161 GeV. 

The three most important uncertainties of this type are: 



The luminosity, which is expected to be known to about 0.5%, including both the known 
theoretical error on the Bhabha cross-section and the expected experimental error, thus 
contributing about 8 MeV to AM^^. 

Unknown higher-order corrections to the theoretical cross-section (see Section 1.1.6 and 
Ref. [jl8[), which if we conservatively assume a value of ±2%, contribute about 34 MeV 



to AMw. 



Uncertainties in the signal efficiency, described in Section 2.3.4 below, which depend on 
the particular decay channel under consideration. 



2.3.2 Background subtraction 

An uncertainty AcTbkgd on the residual background cross-section predicted by Monte Carlo 
propagates as an additive uncertainty to the measured W"*"W~ cross-section, from which the 
background has to be subtracted. It contributes an error 



AM 



w 



da 
dM^ 



w 



AcTbkgd 
eww 



(33) 
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where eww is the signal efficiency, which is found by multiplying the selection efficiency for a 
given channel by the appropriate branching ratio. The quantity \da/dMY^\~^ is shown in Fig. ^ 
as a function of y/s, and is approximately 470 MeV pb~^ at (i/i)"^* — 161 GeV. Experimental 
methods for determining the uncertainty in the background are described in Section |2.3.4| 
below. A systematic contribution to AM-w of about 59 MeV (32 MeV) in the W"'"W~— ^qqqq 
(W+W"— s>qq/i^, W+W^^/z//z/) channels is expected. 



2.3.3 Beam energy and W width 

The error introduced by an uncertainty in the beam energy -Ebeam to Mw is 



AM. 



w 



da 


-1 


da 


dM^ 




dE\]ea,m 



'beam • 



(34) 



In the threshold region the cross-section a^w is essentially a function of the single variable 
^/s — 2Mw only (see Fig. H), and hence the ratio of derivatives in (0) is approximately unity, 
i.e. AMw — 1.0 Aii^bcam- It is estimated that the beam energy will be known to an accuracy 
of 12 MeV. 



The error on Mw introduced by an uncertainty in the W width F-w is 



AM 



w 



da 


-1 


da 


dMy^ 




ciFw 



AFw ^ 0.16 AF 



(35) 



where the value of the ratio of the derivatives corresponds to a/s 
Model, Fw is proportional to the third power of Mw, 



w 



2.080 



M 



w 



80.26 GeV 



161 GeV. In the Standard 
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If the current world average value of Mw = 80.26 ± 0.16 GeV is used then AFw = 12 MeV. In 
contrast, a (combined) measurement of Fw by the CDF and DO collaborations at the Tevatron 
pp collider |^ gives Fw = 2.051±0.061(expt.)±0.013(input) GeV, consistent with the Standard 
Model calculation, but with a much larger error. If we use the Standard Model width, then the 



contribution (2 MeV) to AMw is negligible. See Ref. for a further discussion 



2.3.4 Experimental determination of systematic uncertainties 

Two methods have been proposed to determine the uncertainty in the signal efficiency and 
background cross-section. The ffist method examines the sensitivity of the signal efficiencies 
and background cross-section to uncertainties in fragmentation. Some preliminary studies have 
been performed in the W"'"W~— >qqqq channel by varying the fragmentation parameters Qo, cr^, 
h and Aqcd l^^l within one standard error bounds, and noting the effect on the signal efficiency 



and background cross-section in PYTHIA |^ generated events. 
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The second method uses data and Monte Carlo event samples from LEPl to determine the 
uncertainty in the background cross-section. The selections described in Section |2.2| have been 
scaled to the LEPl centre-of-mass energy and applied to both real and simulated data. The 
difference between data and Monte Carlo gives the error on the background cross-section at 
y/s ~ Mz- The results are then rescaled to ^/s = 161 GeV. It is assumed that the fractional 
uncertainty on the selection efficiencies for the background is the same at the two energies. 

The uncertainty on the signal efficiency appears small from the results of the fragmentation 
test, within the low statistics of the event samples tested. Further studies are needed to quantify 
any effects. Therefore a conservative estimate of the signal efficiency error of 2% is used at 
present, contributing an error of 34 MeV to M^^. The uncertainty on the background cross- 
section in the W+W~— *qqqq channel (8%) is estimated to contribute about 59 MeV to the 
error on M-w- The uncertainty in the W"*"W~^qq/z/ and W"'"W~— >/z//z/ channels, estimated to 
be 50% and 100% respectively, both give an error of about 32 MeV. These errors are expected 
to decrease with further study. 

2.3.5 Conclusions 

Table H summarizes the estimated systematic errors for the W"''W^^qqqq, W+W^— >qq/t/ and 
W^W^ ^Ivlv channels. 



Source 


W+W ^qqqq 


W+W-^qq/z/ 




Luminosity (*) 


8 


8 


8 


HO corrections (*) 


34 


34 


34 


Beam energy (*) 


12 


12 


12 


W width (*) 


2 


2 


2 


Signal efficiency 


34 


34 


34 


Background cross-section 


59 


32 


32 


Total (MeV) 


77 


60 


60 



Table 8: Summary of systematic error contributions, in MeV, to the threshold measurement 
in the W"*"W~— >qqqq, W"'"W~^qq/z/ and W~^W~ ^lulu channels. The quantities denoted (*) 
are common to all channels. The total for each channel is found by combining the sources in 
quadrature. 



2.4 Summary 

Table ^ summarizes the results presented above for the estimated statistical, systematic and 
total errors on Myv (for all decay channels combined) using the threshold method, i.e. by 
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measuring the W"^W~ cross-section at the optimal colhder energy of 161 GeV. Our estimates 
for some of the systematic errors, for example the unknown higher-order theoretical corrections, 
are probably too conservative, and others, for example the uncertainty in the estimates of the 
various background cross-sections, will almost certainly decrease with more study. Nevertheless, 
for the amount of luminosity likely to be available for the threshold measurement the overall 
error is dominated by statistics. 



Total luminosity 


AMw (stat) 


AMw (stat+sysi) 


AMw (total) 


4 X 25 pb"' = 100 pb"' 


134 


139 


144 


4 X 50 pb"^ = 200 pb~^ 


95 


101 


108 


4 X 100 pb~^ = 400 pb"^ 


67 


76 


84 



Table 9: Summary of statistical and systematic errors (in MeV) on Mw from a cross-section 
measurement (all channels) at ^/s = 161 GeV. The total luminosity refers to four experiments 
combined. The third column includes the channel- dependent systematic errors only (see Ta- 
ble and the fourth column includes the overall common systematic error. 



3 Direct Reconstruction of Mw3 

In this section, we discuss the measurement of the W mass by kinematic reconstruction of the 
invariant mass of the W decay products. The statistical precision of this method which could 
be obtained by combining four experiments each with 500 pb~^ at ^/s = 175 GeV, assuming 
100% efficiency and perfect detector resolution, is about 10 MeV, limited by the finite width 
of the W. In practice, this ideal precision will be degraded, partly through loss of statistics, 
but mainly because detector effects limit the resolution on the reconstructed mass. This has 
been studied in detail by the four experiments, using Monte Carlo events with full detector 
simulation. We discuss methods of improving the mass resolution over that obtained by simple 
calculation of invariant masses. In particular, a kinematic fit using the constraints of energy and 
momentum conservation, together with the equality of the two W masses in an event, proves to 
be a very powerful technique for improving the mass resolution, and also turns out to be a useful 
background rejection criterion. For this reason, we concentrate on the channels W^W~^qqqq 
and W^W~— >qq/z/ where the lepton is an electron or muon, for which constrained fits are most 
useful. 

We start by discussing the basic selection of W-pair events in these two channels, and the 
reconstruction of jets. We discuss the techniques of the constrained fit in some detail, followed 
by the determination of Mw from the distribution of reconstructed masses, indicating the 

•^prepared by M. Griinewald, N. J. Kjaer, Z. Kunszt, P. Perez, C. P. Ward 
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statistical error which may be expected. Finally we describe the main sources of systematic 
error pertaining to this measurement. 



3.1 Event selection and jet reconstruction 

The criteria used to select W-pair events are essentially the same as those described in Section 
but at the higher energies used for direct reconstruction the background from Z/7 ^qq is lower, 
so looser cuts can be used. 



3.1.1 W+W-^qqqq 



W"'"W~— >qqqq events are characterised by high multiplicity (about twice that of a Z/7 ^qq 
event at LEPl), high visible energy, and exhibit a four-jet structure. The main background 
comes from [Z^/'j)* — >qq events, for which the cross-section is much higher, as shown in Table |^. 
Typical selection cuts for W"'"W~^qqqq events include: 



high multiplicity of tracks and calorimeter clusters to remove purely leptonic events; a 
significant fraction of W^W^— >qqZz/ and {Z^ /'-/)* — >qq events can also be removed. 

high visible energy and low missing momentum, to remove W~*'W^— s>qqZz/ and radiative 
Z/7 — ^qq events. 

explicit removal of events with an isolated electromagnetic cluster consistent with being 
an initial state photon to remove radiative Z/7 — ^^qq events. 

event shape variables to separate the four-jet W^W~— *qqqq events from the remaining 
background, almost entirely composed of non-radiative Z/7 — >qq events. For example. 
Fig. |11] shows the distribution of y^^, the value of i/cut in the Durham (kx) jet-finding 
scheme at which events change from three jets to four jets, after cuts on the above 
quantities have been applied. The W+W"— s-qqqq events tend to have larger values of 
this variable than the background. Other event shape variables, such as jet broadening 
measures can also be used. 



In Table we show values of efficiency, purity, accepted cross-sections and numbers of 
events produced by typical cuts on these variables. The efficiency of selection cuts tends to 
fall slightly with energy because as y/s is increased the W's are more boosted and event shape 
variables have less separating power. The purity can be further enhanced by using kinematic 
fits, as described below, though at a cost in efficiency. 

In order to reconstruct the W mass, a jet-finder is used to force the selected events to 
contain four jets. Jets are usually reconstructed using both tracks and calorimeter information 
combined to give the best resolution. The typical jet energy resolution is around 20% for 
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Figure 11: Distribution ofy^^, the value ofycut the Durham (k^) jet-finding scheme at which 
events change from three jets to four jets, after cuts on multiplicity, visible energy and missing 
momentum have been applied to events at 175 GeV. The solid histogram shows W+W~— >qqqq 
events, the open histogram background, mainly non-radiative (Z^/j)* — >qq events. The dotted 
line indicates a typical cut value. 
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= 175 GeV 


^ = 192 GeV 




Accepted 


Events 


Accepted 


Events 




cross-section (pb) 


for 500pb-^ 


cross-section (pb) 


for 500pb-^ 


After basic selection cuts: 


W+W ^qqqq 


5.3 


2650 


6.2 


3100 


{TPhY ^qq 


3.1 


1550 


2.1 


1050 




0.1 


50 


0.4 


200 


efficiency 


83% 




79% 




purity 


62% 




62% 




After kinematic 


fit: 








W+W"^qqqq 


4.5 


2250 


5.0 


2500 


(Z77)* ^qq 


1.8 


900 


1.2 


600 


(Z7:0*(z7:0* 


0.08 


10 


0.3 


150 


efficiency 
purity 


71% 
71% 


64% 
77% 


A(Mw) (stat) 


73 MeV 


74 MeV 



Table 10: Typical accepted cross-sections and numbers of events for the signal and main hack- 
grounds for the W^W~— ^qqqq channel for two values of cm. energy, determined from Monte 
Carlo including full detector simulation. Values are given both for basic selection cuts, and 
after demanding a good kinematic fit. The efficiency, purity and expected statistical error on 
for an integrated luminosity of 500 pb"^ are also given. 
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a jet energy of 20 GeV, improving to 15% at an energy of 60 GeV; over this same energy 
range the angular resolution improves from 4° to 1.3° for jets at 90° to the beam direction. 
Studies of various jet finders comparing the reconstructed jets with the initial quark directions 
show no major differences among the commonly used schemes. The W mass may then be 
reconstructed by forming the invariant mass of pairs of jets. For each event, there are three 
possible combinations; Monte Carlo studies show that the combination where the two highest 
energy jets are combined together is rarely the correct one, and combinatorial background can 
be reduced by discarding this combination. The mass resolution can be improved by using 
beam energy constraints or kinematic fits, as described in the next section. 



3.1.2 W+W-^qqlu 

The distinguishing feature of W"'"W~— >qqei/ and W+W"— i>qq/ii/ events is the presence of a 
high momentum, isolated lepton. Typical selection cuts for W+W~— s>qqez/ and W"'"W~-->qq/iz/ 
events include: 

• high multiplicity of tracks and calorimeter clusters to remove purely leptonic events; 
the multiplicity of W~*'W~^qq/z/ events is lower than that of W"*'W~— s>qqqq events, so 
suitable cuts do not remove (Z°/7)* — i>qq background in this case. 

• the presence of a high momentum, isolated, lepton. 

An isolated lepton can be identified as an electron or muon with fairly loose, standard ex- 
perimental cuts with high efficiency (~95%), though only in a hmited acceptance, typically 
I cos(6')| < 0.93. For example electrons can be identified using the match between track mo- 
mentum and energy deposited in the electromagnetic calorimeter, where the pattern of energy 
deposition is consistent with an electromagnetic shower. Muon identification uses matching 
between a track in a central tracking chamber and one in outer muon detectors. The lep- 
ton momentum spectrum and its degree of isolation, as measured by the total scalar sum of 
charged particle momentum plus electromagnetic energy in a 200mrad cone around the lepton, 
are shown in Fig. |12[ Typical efficiencies, purities, accepted cross-sections and numbers of 
events are indicated in Table |ll]. W"'"W~— >qqri/ events can be selected as above, but instead 
of requiring an identified electron or muon, searching for a narrow, low multiplicity jet isolated 
from the other jets. 

The W mass can be estimated from these events by simply forming the invariant mass of 
the hadronic system, scaling to the beam energy, or preferably using the full information in a 
kinematic fit as described below. In this case, the hadronic system is forced to be two jets, which 
are reconstructed as in the W+W~— i>qqqq case. The lepton energy resolution is much better 
than that for a reconstructed jet, as long as care is taken to include all the electromagnetic 
energy associated with an electron. 
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Figure 12: (a) Momentum spectrum of particles identified as leptons and passing isolation cuts 
at 161 GcV. (b) scalar sum of charged particle momentum and electromagnetic calorimeter 
energy in a 200 mrad cone around high momentum identified electrons and muons at 161 GeV. 
(c) and (d) as (a) and (h) for ^/s ~ 175 GeV. In each case the solid histogram shows the 
contribution from leptons in W'*"W~— >qqZi/ events, the open histogram the background. The 
dashed lines show possible cuts. 
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= 175 GeV 




=192 GeV 






Accepted 


Events 


Accepted 




Events 






cross-section (pb) 


for 500pb-^ 


cross-section (pb) 


for 500pb-^ 


After basic selection cuts: 


W+W-->qq/z/ (/ = 


e or /x) 


3.1 




1550 


3.8 




1880 


W+W~^qqri/ 




0.2 




100 


0.4 




200 


(Z77)* ^qq 




0.2 




100 


0.2 




100 


ZO ee 




0.2 




100 


0.3 




150 


efficiency 






77% 






74% 




purity 






83% 






80% 




After kinematic fit: 


W+W-^qq/i/(/ = 


e or /i) 


3.0 




1500 


3.4 




1700 


W+W-^qqri/ 




0.05 




25 


0.06 




30 


(Z77)* ^qq 




0.04 




20 


0.05 




25 


ZO ee 




0.02 




10 


0.05 




25 


efficiency 






73% 






68% 




purity 






96% 






95% 




A(Mw) (stat) 


72 MeV 


93 MeV 



Table 11: Typical accepted cross-sections and numbers of events for the signal and main back- 
grounds for the W^W^^qq/z/ (I = e or fi) channel for two values of cm. energy, determined 
from Monte Carlo including full detector simulation. Values are given both for basic selection 
cuts, and after demanding a good kinematic fit. The efficiency, purity and expected statistical 
error on Mw for an integrated luminosity of 500 pb~^ are also given. 
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3.2 Constrained fit 



After the event has been reconstructed as a number of jets and a number of leptons we next turn 
to the reconstruction of the W mass from these four fermions, which are treated as individual 
objects with measurable quantities. We try to reconstruct the best estimator for the W mass 
on an event by event basis from the measured quantities and the constraints imposed by energy 
and momentum conservation, and the possible additional constraint that the masses of the two 
W's are equal. 

Without imposing any constraints the direct reconstruction of the di-jet mass gives: 



m 



2EiEj{l- cos 9ij), (37) 



where Ei and Ej are the jet energies, 6ij the jet-jet opening angle and where the jet masses 
have been neglected. Taking only the errors on the jet energies into account, as these are much 
larger than the errors on the angle measurements, this then gives: 



a{mij) 



rriij 



\ 




(3J 



Typical jet energy measurement errors of 15% at 45 GeV lead to a relative uncertainty of 10% 
on Mw, and give distributions of reconstructed mass as shown in the top half of Fig. To 
make a precision measurement of Mw it is necessary to improve this resolution by making use 
of the knowledge of the total energy and momentum which is given in an e^e~ collider. 



3.2.1 Rescaling methods 

These methods are especially useful for the analysis of the semi-leptonic channels and most 
have been described previously [||. The basic principle is to write the momentum, energy, 
and equal mass constraints as functions of the measured fermion momenta and solve for those 
variables which have the largest measurement uncertainties, generally jet energies. 

The first step is to include the beam energy constraint, which is equivalent to the constraint 
that the two masses are equal: 

E, + Ej = Eb. (39) 
This leads to a determination of the reconstructed Mw: 

rriij = -^^^pE,E,{l-cose,,). (40) 

Assuming that only the errors on the jet energies are important the error on the mass becomes: 

crjmij) ^ \Ei - Ej 

ITLij Ei + Ej 




(41) 
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Together with the other errors which were neglected, this leads to a relative uncertainty on 
the reconstructed mass of typically 5%. This method is especially suited to the semi-leptonic 
channel in which the leptonic W decay into rvr does not allow advantage to be taken of the 
measured lepton energy. 

In the case of semi-leptonic decays to an electron or muon, we can include the parameters 
of the measured lepton by writing: 

Pu = -EiQi - EiQi - EjQj 
Ei + Ej = Ey, 

E, + E, = Eb, (42) 

where e are unit vectors in the direction of the particles. These five equations with five un- 
knowns, pjy, Ei, and Ej, can be solved, and yield two distinct solutions. This ambiguity leads 
to a problem if the two solutions are close to each other. Taking the solution which is closest 
to the measured jet energies leads to a relative error on Mw of about 4%. 

The two exact solutions of Eq. ( ^21) will give two minima also for the constrained fit in the 
X^(Mw) distribution when the two solutions are close. Monte Carlo studies suggest that about 
40% of the events are afflicted by this problem. Current analyses have not yet included this 
effect in the determination of Mw and one might therefore expect an improvement in resolution 
if this effect is taken correctly into account. 



3.2.2 The constrained fit 

The most effective way to use all the information available in an event is to perform a constrained 
kinematic fit. In such a fit, the measured parameters are varied until a solution is found which 
satisfies the constraints imposed and also minimises the difference between the measured 
and fitted values. Several methods exist to perform this minimization. A traditional one solves 
the problem using Lagrange multipliers, minimising 

S{y, A) = {y-y-yY-\y-yo) + 2A • /(y), (43) 

where V is the error matrix, y the fitted variables, y^ the measured values, A Lagrange multi- 
pliers, and f{y) the constraints written as functions which must vanish. An equivalent method 
is to use penalty functions where terms of the type jcP' are added to the x^- The procedure 
then minimizes the total in an iterative way, for each step decreasing the o of the penalties. 
Results are in general the same but the convergence is slightly slower. 

The inputs to the fit are measurements of the energy and angles, of the four jets in the 
hadronic final state or of the two jets and lepton in the semi-leptonic final state, together with 
their error matrix. Errors on jet parameters can be extracted from data or Monte Carlo, and 
may be functions of both energy and position measured in the detector. In practice most of 
the jet measurement errors are nearly uncorrelated, making the error matrix diagonal. For the 
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jet masses two different strategies are used. Either the jets are assumed to be massless and the 
measured jet energy is used as the measured jet momentum, or one includes the reconstructed 
jet masses in the fit. (Note that the r lepton can be treated in nearly the same way as jets, 
except that the errors on the transverse momenta are given by the missing mass and that the 
r mass is used explicitly in the fit. This has not yet been studied in detail.) 

The fit can be performed using only the constraints of energy and momentum conservation, 
or also including the constraint that the masses of the two reconstructed W's be equal. In 
the hadronic channel, this gives a 4C or 5C fit respectively. In the semi-leptonic channel, 
the number of constraints is reduced by three because the parameters of the neutrino are 
unmeasured, resulting in a IC or 2C fit. The equal mass constraint can be included exactly, or 
the width of the W can be taken into account by adding a term to the proportional to the 
difference in mass of the two W's. In practice, because the mass resolution is larger than the W 
width, both methods give almost the same results. If an equal mass constraint is not applied, 
the reconstructed masses of the two di-fermion systems are strongly anticorrelated. Thus only 
the average invariant mass can usefully be extracted per event. The inclusion of an equal mass 
constraint is preferred over a fit using only energy and momentum conservation because it gives 
improved mass resolution and superior background rejection. 

In the 4-jet channel we do not know which jets are to be combined to produce the heavy 
particles we are interested in reconstructing. Therefore the constrained fit is performed for 
all three possible combinations. In the case with no equal mass constraint there are three 
different mass solutions with the same x^- When an equal mass constraint is imposed the three 
different solutions will in general have different and one can therefore use this information to 
distinguish between the solutions. The procedure employed by most analyses is to first define 
a mass window, wide enough not to bias the mass measurement, and if more than one solution 
exists inside this window choose the solution with the lowest x^- This procedure is however not 
perfect and one is left with a fraction of wrong combinations as shown in Fig. |T3|. 



3.2.3 Results of the fit 



In Fig. |l3|we show invariant mass distributions before and after the fit, in the latter case taking 
only those events which give a good fit. Before the fit, the mass distributions are very broad. 
After the fit the mass resolutions obtained are typically 3.5% for the W"''W~— >qqZz/ channel 
and around 2.5% for the 4-jet channel. It is also clear from this figure that the of the 
constrained fit can be used to eliminate possible background which does not comply with the 
W'^W" hypothesis: the level of background is much lower in both channels after demanding a 
good fit. However, a fraction of W"'"W~ events, varying from about 10% to 30% depending on 
the channel, also fails to give a good fit. These events are discussed below. Typical values of 
efficiency and purity after the fit are given in Tables and 

From the constrained fit we can calculate the error on the fitted mass. This error is highly 
correlated to the actual value of the fitted mass, so selecting events with a simple cut on this 
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Figure 13: Distributions of the invariant mass for selected events for the W"'"W ^qqqq channel 
at ^ = 192 GcV (left) and the W+W"^qq/z/ channel at ^ = 175 GeV (right). The 
top plots are before the kinematic fit, the bottom ones after a five- constraint fW+W~— >qqqq^ 
or two- constraint (W'^W"— >qq/z/^ kinematic fit. For the W"^W~^qqqq channel, the open 
histograms show the correctly found jet-jet combinations, the singly hatched areas correspond to 
the incorrectly found jet-jet combinations and the cross-hatched parts indicate the background 
from other sources (mainly qq{'^)). For the W"''W^^qq/z/ channel the open histograms show 
signal events, the cross-hatched areas background including the W"'"W~— >-qqTi/ contribution. 
The statistics correspond to an integrated luminosity of 500 ph~^ . 
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quantity would seriously bias the measurement. The reason for this effect is related to the 
kinematics of the events. When the mass approaches the kinematic limit the precision on the 
sum of the unknown masses will be better and better while the precision on the difference be- 
tween the masses will deteriorate. This leads to a mass resolution that broadens with increasing 
energy. Going from 175 GeV to 192 GeV typically increases the measurement errors by 25% 
on an event by event basis. The effect of this on the final statistical error on Mw is diluted by 
the fixed Fw and compensated by an increase of the W''"W~ cross-section. 

The constrained fit assumes that the errors on the measured quantities arc Gaussian and 
uncorrelated between jets. Several effects lead to non-Gaussian errors and correlations. Each of 
these will lead to tails in the distributions and hence to a peak at low probability for the fitted 
X^. The most important is gluon radiation, but also overlapping jets, initial state radiation, Fw, 
and acceptance effects play a role. The hard gluon radiation is of course in direct disagreement 
with the treatment of jets as independent objects. Even rather soft gluon radiation lead to jets 
being broadened in a specific direction, giving correlations that are not included in the current 
implementations of the constrained fit. Studies have been performed to try to recover some of 
the 4-jet events which fail to give a good fit by treating them as 5-jet events. However, although 
some fraction of these events then give a good fit, the combinatorial problem is severe, and it 
appears that their inclusion has little effect on the ultimate mass resolution. 

3.2.4 Inclusion of initial state radiation 

As will be seen below, energy lost in initial state radiation biases the fitted mass if it is not 
included in the fit. Initial state radiation can be included in the constrained fit using the 
following procedure. We know that there is a large probability that the photon is produced 
close to the collision axis and hence not detected. We can therefore as a good assumption take 
the momentum to be collinear with the 2;-axis. We also know to a high precision the expected 
distribution of the photon energy. In a simple approximation this reads: 



where x — E^j ^/s and /3 is a parameter which is smaller than unity and depends on ^/s. Since 
this is non-Gaussian we cannot take this term directly into the expression. Instead we 
introduce the likelihood concept and use the standard -2 In(likelihood) as the term to add to 
the x^. If we just use the probability as the likelihood this approach will not work since the 
distribution has a pole for a; — > 0. Instead one can choose to use the confidence limit as an 
estimator for the likelihood: _ 



With this formulation the constrained fit works but its implementation is rather difficult since 
one has to divide the fit into two parts: one where one assumes the photon goes in the forward 

direction and one where one assumes the opposite. When the fitted approaches zero the 
first order derivative of will still diverge, but this only means that the fit prefers the zero 
solution, since the measurement does not have sufficient resolution to distinguish between no 



p{x)dx — (3x^ ^dx, 



(44) 




(45) 
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7isR and a 7isr with a small energy. Monte Carlo studies show that photons below about 3 GeV 
cannot be resolved and when photons are above typically 8 GeV the W bosons can not be on 
mass-shell. The final improvement in the Mw measurement is therefore limited. 



3.3 Determination of the mass and width of the W 

In this section we describe several strategies for extracting Mw from distributions of recon- 
structed invariant masses such as those show in Fig. As we aim for a precision measurement 
of Mw with sub-permille accuracy, this is a non-trivial task, because we have to control any 
systematic effect to an accuracy of a few MeV. The total width of the W boson, Fw, may ei- 
ther be extracted simultaneously with the mass, or the functional dependence Fw = Fw(Mw) 
of the Standard Model may be imposed as a constraint for increased accuracy on Mw- The 
methods to analyse the data in terms of Mw and Fw fall into four groups: (1) Monte Carlo 
calibration of simple function; (2) (De-) Convolution of underlying physics function; (3) Monte 
Carlo interpolation; (4) Reweighting of Monte Carlo events. 



In general, all methods make use of Monte Carlo event generators [|T^] and detector simu- 
lation to determine the effects of the detector such as resolution. Thus any method has to be 
checked for possible systematic biases introduced by using Monte Carlo event samples gener- 
ated with certain input values for Mw and Fw- In addition, systematic errors may arise due 
to deficiencies in the Monte Carlo simulation describing the detector and/or the data. Other 
systematic errors arise from the technical implementation of the fitting methods, such as fit 
range, bin width in case of binned data, choice of functions to describe signal and background 
etc. 



3.3.1 Monte Carlo calibration 

To fit the invariant mass distribution, this method uses a simple, ad-hoc function, e.g., a double 
Gaussian or a Breit-Wigner convoluted with a Gaussian, to describe the signal peak and another 
simple function to describe the background. One of the fit parameters is used as an estimator, 
M, for the W mass, e.g., the mean of the Gaussian or the Breit-Wigner. The same fitting 
procedure is applied to both data and Monte Carlo events. Since for the latter the input W 
mass, M^*^, is known, the Monte Carlo result is used to evaluate the bias A of this method, 
A = M^^ - M^^, where this bias may depend on the final state analysed. The mass of the W 
measured in the data is now simply given by the estimator M'^'^^^ derived from fitting the data 
distribution, corrected for the bias A evaluated from Monte Carlo events: Mw = M'^^^^ — A. 

This procedure automatically takes into account all corrections for all biases as long as they 
are implemented in the Monte Carlo simulation, such as initial-state radiation, background 
contributions, detector resolutions and efficiencies, selection cuts etc. The knowledge of how 
well the Monte Carlo describes the underlying physics and the detector response enters in the 
systematic error on the bias A. 
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The fundamental drawback of this method is that the simple function used in the fit is not 
unique. Depending on its choice, even different statistical errors on Mw can be obtained. In 
addition, the estimate of the bias correction A depends itself to some extent on the Monte 
Carlo parameters M-^*^ and T^^, or on the centre-of-mass energy, ^/s, A = A{s; M^'^,T^^). 
Such a dependence, however, can be corrected for by iteration. 



3.3.2 Convolution 

The drawbacks listed above are alleviated in the convolution method. Here, the correct func- 
tion, i.e., the underlying physics function, is used as a fitting function. This function is simply 
the differential cross-section in the two invariant masses (denoted by mi, 7712). Note that this 
function is not a simple Breit-Wigner distribution in mi and m2 due to phase space effects 
and radiative corrections. Several analytical codes exist (e.g., GENTLE etc. ||18|), which calcu- 



late the differential cross-section including QED corrections, '^('^'^^w.rw) function of the 

° ^ ' dmiam2 ' 

centre-of-mass energy, y/s, and the Breit-Wigner mass and total width of the W boson, Mw 
and Fw- 

The effects of the detector are included by convolution. The prediction for the distribution 
of the reconstructed invariant masses (denoted by mi,m2) is thus given by: 

dV(s;Mw,rw) [ A [ A ni ^ dV(g;Mw,rw) 

— — — = I ami I 6.1712 G(s\ nil, 1712, nil, 1^2) ■ -, . (46) 

ami am2 J J ami am2 

The transfer or Green's function G, which depends on the final state analysed, can be inter- 
preted as the probability of reconstructing the pair of invariant masses (mi, given the event 
contained the pair of true invariant masses (mi,m2). Several simplifications for G are possible, 
down to a 1-dimensional resolution function G = G ([mi + m2]/2 — [mi + m2]/2). 

The actual fitting of invariant mass distributions can be performed either on the measured 
or on unfolded distributions. In the former case, the underlying physics function, d^m^dm^^^ 
is convoluted with G, and the result is fitted to the data with and Fw as fit parameters. 
In the latter case, the measured distribution, '^(M^w.rw) jg ^^.^^ unfolded for detector effects, 

' ' dmidm2 ' ' 

by applying the "inverse" of G. The underlying physics function, '^dnndrr^^^ ' '^^^ ^'^^ 
fitted directly to the unfolded distribution. From a purely mathematical point of view, both 
methods are equivalent. In practice, however, G is determined only up to a certain statistical 
accuracy. Since folding is an intrinsically stable procedure in contrast to unfolding, which is 
unstable, the former method is preferred. Reference ||4^ gives more details on the features of 



unfolding procedures. Backgrounds are described by simple functions, derived from data and 
from Monte Carlo simulations, which are added to the signal. 

This method allows cuts on generated invariant masses (mi,m2) and ISR energy loss, since 
these are the only variables used in most semianalytical codes. In addition, cuts on the re- 
constructed invariant masses (mi,m2) are possible. Since cuts on other variables cannot be 
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applied, it must be checked whether selection cuts bias the invariant mass distribution, for 
example by testing the method with Monte Carlo events (cf. method (1)). 

3.3.3 Other Monte Carlo based methods 

The problems of the previous two methods can be solved by a Monte Carlo interpolation 
technique. Several samples of Monte Carlo events corresponding to different input values of 
M^c and r^*^ are simulated, e.g., in a grid around the current central values for (Mw,rw) 
extending a few times the total (expected) error in all directions. The Monte Carlo samples 
of the accepted events are compared to the accepted data events, thereby taking the influence 
of event selection cuts into account. Backgrounds are included by adding the corresponding 
Monte Carlo events. The compatibility of the invariant mass distributions is calculated, e.g., 
in terms of a quantity. Interpolation of the within the generated (Mw^, T^^) grid allows 
to find the values Mw and Fw which minimise the y^. Like method (1), this method corrects 
automatically all possible biases due to all effects considered in the Monte Carlo simulation. 
The only drawback is that a rather large amount of Monte Carlo events must be simulated. 

This problem, however, can be solved by a reweighting procedure. In that case only one 
sample of Monte Carlo events is needed, which has been generated with fixed values M^'~^ 
and r^^. Event-by-event reweighting in the generated invariant masses (mi, 7712) is performed 
to construct the prediction for the invariant mass distributions corresponding to values 
and r^. These distributions are then fitted to the data distributions. Since individual Monte 
Carlo events are reweighted, it is straight forward to implement the effects of selection cuts. 
Moreover, using a Monte Carlo generator also in the calculation of the event weights, it is even 
possible to extend further the set of variables on which the event weights depend to include any 
kinematic variable describing the four-fermion final state, such as the reconstructed fermion 
energies and angles. 

3.3.4 Expected statistical error on Mw 

So far the experiments have studied methods (1), (2) and (3) |^ 0. For the statistical 

errors on Mw quoted in Tables and 0, the experiments have mainly used method (1), which 
is adequate for this purpose. It should be noted that the more involved analyses (2), (3) and 
(4) do not aim for a reduction in the statistical error on Mw. This error is fixed by the number 
of selected events, the natural width of the W boson and the detector resolution in invariant 
masses. Instead the advantage of these methods lies in the fact that they allow a thorough 
investigation of the various systematic biases arising in the determination of the W mass. It is 
possible to disentangle systematic effects due to background diagrams, initial-state radiation, 
event selection, detector calibration and resolution etc. Thus, in order to have a better control 
of systematic effects, it is envisaged that the final analyses will use the more involved strategies 
(2), (3), (4) or a combination thereof. 
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As shown in Table ITDI, the statistical error on Mw from the W^W~— i>qqqq channel is 
expected to be about 73 MeV, roughly independent of y/s, for an integrated luminosity of 
500 pb~^. For the W'^W^— s>qq/i/ channel a similar value is expected at 175 GeV, but at higher 
energies the worsening resolution on Mw causes this value to increase somewhat. 



3.4 Systematic errors 

Three classes of systematic error have been envisaged: errors coming from the accelerator, from 
the knowledge of the underlying physics phenomena, and from the detector. The last two classes 
are sometimes related as some physics parameters may affect the detector response. In contrast 
with the threshold cross-section method, effects which distort the mass distribution must be 
considered here. Most of the following error estimates have been obtained using simulations. 
The models used will be checked against LEP2 data when available and represent today's state 
of the art. 



3.4.1 Error from the LEP beam 

The direct mass reconstruction method relies on constraining or rescaling the energies of the 
reconstructed W's to the beam energy. The error on the beam energy is foreseen to be less than 



12 MeV (see Section |1 . 1| ) , translating to 12 MeV x M^^ / E},f,am on the mass. The beam energy 
spread will be ~ 40 MeV x s/(91 GeV)^ 0] with a negligible effect on the statistical error. The 
operation with bunch trains should not significantly affect the beam energy. Experience gained 
at LEPl will guarantee a very good follow up of the behaviour of the beam even at short time 
scales. 



3.4.2 Error from the theoretical description 

The most important distortion of the mass distribution comes from initial state radiation (ISR). 
The average energy carried by radiated ISR photons is {Ejsr) = 1.2 GeV per event at i/i = 
175 GeV and rises linearly up to 2.2 GeV at ^/s = 200 GeV. From the beam energy constraint, 
one expects an average mass shift of the order of {Eisr)M-w/ y/s, i.e. about 500 MeV at 
175 GeV. A fit to the mass distribution gives more weight to the peak yielding a shift of only 
200 MeV at this energy. As a direct measurement of the ISR spectrum seems impossible to 
an accuracy relevant for the W mass measurement, we will rely on the theoretical calculations 
as described in Section |1.6.5| which present an error of 10 MeV on {Ejsr). Together with 



uncertainties on the shape of the spectrum, this translates into an error smaller than 10 MeV 
on the W mass. 

Background events contained in the final sample will also distort the mass distribution 
(Fig. |T^). These backgrounds can be monitored with the data themselves under a dedicated 
analysis. The most important source of background comes from (Z^/j)* — >qq events which has 
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a flat mass distribution in the fit area. The effects of an error on the level or in the shape of 
this background have been estimated by changing its production cross-section by 20% and by 
shifting the background mass distribution by 0.5 GeV resulting in an error of 12 MeV for the 
W"''W~^qqqq channel and 6 MeV for the W"*"W~— *qq/z/ channel. The background four jet 



rate as measured at LEPl is not perfectly reproduced by the simulations |52| and will have to 
be better determined. 

The fragmentation model for W events is the same as was fitted to Z events at LEPl. To 
estimate errors from this model, each of its parameters has been varied and the error is derived 
according to a 1 standard deviation from the best fit value. The total error is 16 MeV, which 
should be taken as an upper limit as the parameters have been varied separately without taking 
into account correlations. 

Colour reconnection and Bose-Einstein correlations may seriously affect the mass distribu- 
tion in the W"''W~^qqqq channel. These are discussed in detail in Section H. 



3.4.3 Error from the detector 

Effects of miscalibration appear to be very small. This is due to the use of the beam energy 
which gives a precise scale. For instance, a conservative shift of 2% in the momentum of charged 
tracks, and miscalibration of 5% for electromagnetic and hadronic calorimeters produces a shift 
of 10 MeV on the W mass. 

During the LEP2 era, it is foreseen to run for short periods at -^i = Mz in order to calibrate 
and determine the efficiencies of the detectors with high statistics. Genuine LEP2 events like 
e"'"e~ — *• Z°7 will also be used to monitor the hadronic recoil mass to the photon. Two Z events 
produced at LEPl can be mixed and boosted to a LEP2 energy to study the adequacy of the 
mass reconstruction in comparison with the simulations. 

In the W"''W~-^qqqq channel, the wrong assignment of a jet to a W distorts the mass 
distribution giving rise to a low mass tail. The amount of such misassignments depends on 
the algorithms used. By changing parameters in such algorithms we estimate an error of up to 
5 MeV. 

Several methods to extract the W mass from the measurements have been studied (see 



Section |3.3| ). If a simple function is fitted to the measured mass distribution, errors from the 
choice of the mass fit range are to be expected, and are typically 10 MeV. A convolution 
of the physics function with (or a deconvolution from) detector effects should largely reduce 
such errors. However this convolution depends on the accuracy of the simulation of the event 
structure, for instance on the fragmentation model used (see above). 

In the present simulation studies, limited statistics have been generated. The statistical 
error on the mass shift with 500k events simulated is of the order of 10 MeV. 
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12 


6 


calibration 


10 


10 


MC statistics 


10 


10 


mass fit 


10 


10 


jet assignment 


5 




interconnection 


? 




total 


31 


29 



Table 12: Estimated systematic errors on Mw J'er experiment in MeV. Colour reconnection 
and Bose-Einstein effects are not included in the total. 



3.4.4 Common errors 



The effects mentioned above are summarized in Table [T^ and total to about 30 MeV for each 
channel at =175 GeV. Although errors from the beam energy, ISR and background will 
vary with energy, we do not expect large variations over the 165 to 190 GeV range. 

Errors on the beam energy and ISR are common to all experiments. Errors from fragmen- 
tation and background are partially common as far as the physical parameters are concerned, 
but a particular analysis will retain a level of background or be sensitive to fragmentation tails 
in a different way from another analysis, hence with a slightly different error. If these errors 
were considered common to the four experiments, they would amount to a total of 25 MeV 
(23 MeV) in the W+W~^qqqq (W+W-^qq/z/) channel. 



3.5 Summary 



In Table |T3| we show the total error on Mw which might be expected by combining four 
experiments each with an integrated luminosity of 500 pb~^ at y/s = 175 GeV. As discussed 
in the previous section, systematic errors from the beam energy measurement and initial state 
radiation, fragmentation and backgrounds are considered common to both channels and to all 
experiments. The total expected error is around 34 MeV. with roughly equal contributions 
from statistics and systematics. This value is expected to rise, but only very slowly, with i/i 
because of the worsening mass resolution. 
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Source 


W+W"^qqqq 


W+W-^qq/z/ 


Combined 


Statistical 


36 


36 


25 


Common systematic 


25 


23 


23 


Uncorr. systematic 


9 


9 


6 


total 


45 


44 


34 



Table 13: Estimated total error on Mw, in MeV, which could be obtained by combining four 
experiments each with an integrated luminosity of 500 pb~^ at y/s = 175 GeV. Colour recon- 
nection and Bose-Einstein effects are not included in the systematic error. 

4 Interconnection EfFects|^ 



4.1 Introduction 

The success of the precision measurements of the W boson mass Mw strongly relies on accurate 
theoretical knowledge of the dynamics of the production and decay stages in e"'"e~ — > W'^'W^ 
4 fermions. Owing to the large W width, Fw, these stages are not independent but may 
be interconnected by QCD (and electroweak) interference effects, which must be kept under 
theoretical control. Interconnection phenomena may obscure the separate identities of the two 
W bosons, so that the final state may no longer be considered as a superposition of two separate 
W decays [|5T], Thus the "direct reconstruction" method of measuring Mw at LEP2 using 
the hadronic (qiq2 q3q4) channel has an important caveat — the colour reconnection effects 
Q — which may distort the mass determination |Q. Another delicate question is 
whether Bose-Einstein effects could induce a further uncertainty in the mass determination 



pi p 



The hadronic effects mentioned above are all well-known from other processes. Interferences 
in the production and decay of unstable objects are observed e.g. in multipion final states at 
low invariant masses QCD interconnection could be exemplified by J/^ production in B 
decay: in the weak decay b cW^ ^ ccs the cc J/ip formation requires a "cross-talk" 
between the c-|-s and the c-|-spectator colour singlets. Bose-Einstein effects are readily visible in 
high-energy multiparticle production, including LEPl. In view of the precedents, the working 
hypothesis must be that interconnection effects are at play also for W"'"W~ events. 

The space-time picture of hadroproduction in hadronic W'''W~ decays plays a very impor- 



tant role in understanding the physics of QCD interference phenomena |]52| , |56| . Consider, for 
instance, a centre of mass energy of 175 GeV where the typical separation of the two decay 
vertices in space and in time is of the order of 0.1 fm. A gluon with an energy u ^ Fw therefore 
has a wavelength much smaller than the separation between the W"*" and W~ decay vertices, 
and is emitted almost incoherently either by the qiq2 system from one W or by the q3q4 one 

* prepared by V.A. Khoze, L. Lonnblad, R. M0ller, T. Sjostrand, S. Todorova and N.K. Watson. 
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from the other. Only fairly soft gluons, uo ^ Tw, feel the joint action of all four quark colour 
charges. On the other hand, the typical distance scale of hadronization is about 1 fm, i.e. much 
larger than the decay vertex separation. Similarly, observed Bose-Einstein radii are of this 
order. As a result, the hadronization phase may induce sizeable interference effects. 

As one could anticipate from the above, the perturbative QCD interconnection effects appear 
to be very small In examining the size of the non-perturbative W mass distortions one 
has to rely on existing QCD Monte Carlo models rather than on exact calculations. These 
models have done a very good job in describing a large part of the experimental data at the 
Z" resonance, so it seems plausible that they (after the appropriate extensions) could provide a 
reliable estimate of the magnitude of the interconnection-induced shift in the W mass. However, 
one has to bear in mind that there is a true limit to our current understanding of the physics 
of hadronization |5^, ^ , so the actual value of the shift may be beyond our current reach. 



One of the achievements of the activity of our working group is that the colour reconnection 
physics has been tested in several approaches [^, ^ |5^. The Bose-Einstein studies 



are somewhat lagging behind, but progress is visible also here |]5^, |60|, 



There is another challenging reason to study the phenomenon of colour recoupling in 
hadronic W'''W~ events: it could provide a new laboratory for probing non-perturbative QCD 
dynamics |^ . The very fact that different assumptions about the confinement forces may give 
different predictions for various final-state characteristics means that it might be possible to 
learn from experiment about the structure of the QCD vacuum l5l|, |5^, |53|. 



Finally, note that there are other effects — this time originating in purely QED radiative 
phenomena — which, in principle, prevent the final state from being treated as two separate 
W decays. For instance, final-state QED interactions, for the threshold region exemplified by 
the Coulomb forces between two unstable W bosons, induce non-factorizable corrections to the 
final-state mass distributions in e"'"e~ W"*" W~ — >■ 4 fermions |]32| , |33| , |3^ . Of course, there is 
no reason why all such QED effects cannot, in principle, be computed to arbitrary accuracy, and 
taken into account in the mass determination. However, at the moment no complete formulae 
are available which could be relevant for the whole LEP2 energy range. This topic certainly 
needs further detailed studies and a comprehensive activity here is in progress. For further 
comments see Section 0. 



4.2 Colour reconnection 

Several quite different colour reconnection models have been presented by now. A crude survey 



of features is found in Table [13. The listing is by no means complete; several simpler toy models 



have also been proposed The three main philosophies are the following: 



A colour-confinement string is created in each W decay, spanned between the decay- 
product partons. These strings expand out from the respective decay vertex and even- 
tually fragment to hadrons. Before that time, a reconnection can occur by a space-time 
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encounter between sections of the two strings. In one extreme (QCD vacuum structure 
analogous to a type I superconductor) the reconnection probabihty is related to the in- 
tegrated space-time overlap of the extended strings, in another (likewise, with type II 
superconductor) to the crossing of two string cores (vortex lines). Many subvariants are 



conceivable. The models in |59| represent further refinements of the ones presented in 



Perturbative QCD prefers configurations with minimal string length in TP decays. Here 
length is defined in terms of the A measure, which may be viewed as the rapidity range 
along the string: A ~ I]i lii("^?/''^p); where rrti is the invariant mass of the i'th string 
piece and rup sets a typical hadronic mass scale. It is plausible that, when the partons of a 
W pair are separating and strings formed between them, such configurations are favoured 
which correspond to a reduced total string length. Therefore a reasonable criterion is to 
allow A-reducing reconnections at the scale of Fw, within the limits of what is given by 



colour algebra factors. The models in |^ and |^ are variants on this theme. 



• In a cluster approach to fragmentation, clusters are formed by the recombination of a 
quark from one gluon branching with an antiquark from an "adjacent" branching (alter- 
natively primary quark flavours or diquarks). Normally "adjacent" is defined in terms of 
the shower history, but another reasonable criterion is the space-time size of the formed 
clusters. Therefore reconnection relative to the ordinary picture could be allowed anytime 
the sum of the squared sizes of the formed clusters is reduced. Restrictions are imposed, 
so that e.g. the quark and antiquark from a gluon splitting cannot form a cluster together. 
This approach is found in . 

The models need not be viewed as mutually contradictory. Rather, each may represent some 
aspect of the true nature of the process. 



The main philosophies can be varied in a multitude of ways — Table |14| is to be viewed 
as an appetizer to a proper appreciation of this. Program details can be found in the QCD 
event generators section and in the original literature. The overall reconnection probability 
depends on free parameters in all the models, since we do not understand the non-perturbative 
dynamics; colour suppression factors like l/A*"^ may be present, as in the perturbative phase, but 
are multiplied by unknown functions of the space-time and momentum variables. Furthermore, 
the models can only be incomplete representations of the physics, e.g., the effects of "negative 



antennae/dipoles" |£2| and of virtual-gluon corrections remain to be addressed. 

Several studies have been performed to estimate the effect (bias) of colour reconnections 
on the W mass measurement, and to evaluate the possibility of determining experimentally 
whether such phenomena are taking place in the W^W~ system. Reconnection effects could 
also give visible signals at LEPl. To allow unbiased comparisons, each reconnection model 
should be retuned to general event-shape data with the same care as for the "no reconnection" 
scenarios. This process is now under way. The Ariadne studies make use of a retuning of the 
model with reconnection effects, by its author [0, |6^ . 
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authors 

I cicl cIlCc 


Khoze 
Sjostrand 

m 


To dor ova 
59 


Gustafson 
Hakkinen 

[53, 36 


Lonnblad 

|57| 


Webber 
If; si 


based on 


Pythia 


Ariadne 


HERWIG 


I ecoiiiiectioii 
criterion 


space-time overlap (I) 
or crossing (II) of strings 


string length 
reduced 


Cluster space— tune 
sizes reduced 


reconnection 
probability 


I: free parameter 
II: partly predicted 


free 
parameter 


partly 
predicted 


free 
parameter 


model of 
all events 


yes 


yes 


no 


yes 


yes 


space-time pictu 

W vertices 
parton shower 
fragmentation 


re implemei 

yes 
no 
yes 


ited for 

yes 
yes 
yes 


no 


(- 

no 


- = not applicable) 
yes 
yes 


multiple 
reconnections 


no 


yes 


yes 


yes 


yes 


reconnection 
inside W/Z 


no 


yes 


yes 


yes 


yes 


change of event 
properties 


a 

in^ 


most 
i^isible 


small but 
visible 


visible, needs 
retuning 


large, needs 
retuning 



Table 14: Survey of colour reconnection models. The information should be taken as indicative; 
only the original literature gives the ideological motivation behind the choices. 
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In order to set limits on possible W mass shifts, the experimental sensitivity with which 
reconnection may be observed has to be evaluated for each model (as a function of the model 
parameters, where appropriate). If no effects are observed in the LEP2 data, this sets a limit 
on the reconnection probability, which can be turned into an estimate of the systematic uncer- 
tainty in the W mass. To complicate matters, different models will in general predict different 
mass shifts for the same probability of reconnection. This complication will also appear if 
reconnections are observed in the data and one would like to correct the observed W mass 
accordingly. 

Predictions for the systematic error in W mass measurement are presented in Table 
It is important to understand that the value of the W mass shift depends very strongly on 
the definition of this shift, as well as on the method used to reconstruct the jets, the fitting 
function used, and the tuning of the models. These factors need to be controlled strictly to 
allow meaningful comparisons. Three estimates of the mass shift have been studied: 



• Averaging Some Monte Carlo authors |^2|, ^ form a distribution (difference between 
reconstructed and generated W mass) event-by-event, and use the difference in the means 
of two such distributions (with and without reconnection) as an estimator of the system- 
atic shift. 

• Fitting Two mass distributions are formed, one with and one without reconnection, 
and the W mass is obtained by fitting each separately. The function used was a Breit- 
Wigner, modulated by phase space factors. The difference between the fit results gives 
the estimated mass shift. The definition of jets and correspondence between jets and W's 
followed method 3 of 

• Detector The analysis was performed as if using real data. Mass distributions are 
obtained after the events (with initial state radiation) have been passed through a LEP 
detector simulation program, thereby including typical experimental acceptance and qual- 
ity cuts. The W mass is determined for each event using a 5 constraint kinematic fit (see, 
e.g.. Section H). The shift is taken as the difference between two fitted mass distributions 
as above, although the definition of jets, the association between jets and W's and the 
fitting function were different. The experimental results presented here were obtained 
using a simulation of the OPAL experiment [|63| . 



Table |15] contains an (incomplete) comparison of the shifts observed using all three definitions 
of the shift in the W mass, where the errors are from finite Monte Carlo statistics. As indicated, 
the three are not expected to agree. It is noted that there is an upwards shift after detector 
simulation. The study of including the same kinematic fitting, jet and W reconstruction, 
was repeated using events without detector simulation. Mass shifts which are l-2a (statistical) 
lower than the detector simulation fits were observed; they were stable against changes in 
fit range, and cuts on Pt > 150 MeV/c and kinematic fit probability (>1%). They are also 
consistent with the shifts in the 'fitting' column in Table |T5|. 
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Model 


Reconn. 
prob.(%) 


AMw (MeV) 


Sensitivity (cr) 


averaging 


fitting 


detector 


Central 


Interjet 


|5| Type I, p=0.9 


46.7 


27±14 


61±20 


130±40 


2.4 


1.9 


p=0.6 


36.5 


22±14 


68±20 


80±40 


2.0 


1.5 


p=0.3 


22.4 


10±14 


35±20 


60±40 


1.4 


1.2 


Type II' 


34.2 


-22±14 


-25±19 


50±40 


0.8 


0.8 


H Type I, p=l. 


22 


0±14 


33±18 


- 


- 


- 


Type I, p=100. 


64 


17±14 


47±18 








Type II, (i=0. fm 


8 


-18±14 


-1±18 








Type II, rf=0.01 fm 


34 


-2±14 


32±18 








11 


100 


70±20 




150±40 


6.5 


6.8 


m 


52 


30±10 


58 ± 21 


154±36 






13 retuned(|6|]) 


46 


-6±14 


17 ± 19 








158| 


8 




11±10 









Table 15: Summary of effects on Mw at 175 GeV. See the text for details. The p parameter 
in Type I models relates the string overlap integral to the reconnection probability. The d 
parameter gives the vortex line core diameter in Type II models. The numbers for can be 
rescaled by the reconnection probability, which here is a free parameter. The sensitivities in 
the last two columns include detector simulation [|63i. Several of the numbers in the table are 



preliminary and further studies are underway to understand them. 



If the shift is genuinely caused by specific experimental cuts, it may be possible to redesign 
those cuts to reduce AMw to the level seen without detector simulation. If so, we might also 
construct a related observable for the detection of reconnection. 

Two related experimental methods of detecting the presence of reconnection phenomena at 
LEP2 have been studied, including OPAL simulation 
bin and interjet multiplicity . 



63]: multiplicity in a central rapidity 



In ref. a significant decrease was found in the central rapidity region for events with a 
thrust cut T > 0.76. (A larger difference but poorer statistical significance was obtained for 
more restrictive thrust cuts.) An independent, PvTHlA-based implementation of the GH model 
[p3| showed a much reduced signal. It was here pointed out that one reason for the difference is 
the effect of W polarization, which was neglected in ref. [^. Recent studies confirm that 
the polarization is important; it reduces the signal by roughly a factor of 2 compared to the 
original results. However, although significantly reduced, a clear signal for reconnection effects 
is still visible in this appears to be larger than the result obtained in ref. 



Following the prescription of , the central multiplicity is defined as the number of parti- 
cles, n, observed within a given rapidity interval. In order to quantify the level at which colour 
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reconnection could be observed, the sensitivity at a given multiplicity, n, is defined as [p3| 



^ ^ -^rec 



1 



/ 



1 + Xl^norm(0 



(47) 



where A^rec(0 ^-nd A^norm(^) are the number of events in the reconnected and normal event sam- 
ples having a central multiplicity i, respectively. The sensitivities given in Table |15| are obtained 
assuming that 3 000 four quark events have been observed, that is about what is expected from 
the nominal luminosity of 500 pb~^, by scaling the sensitivities from samples of 50 000 generated 



events by y3 000/50 000. If the results of the four experiments are combined, the statistical 
significance will double. With such statistics a tentative conclusion from the models so far 
studied is that reconnection effects could be observed down to 25-30%, corresponding to a 3a 
effect. A smaller integrated luminosity is not very encouraging at an energy of 175 GeV. It 
should also be remembered that the significance level is based on an assumed perfect knowledge 
on the no-reconnection multiplicity distribution; an entirely experimental procedure based on 
a comparison with mixed leptonic-hadronic events would have a smaller sensitivity. 

A reconnection probability of 30% could induce a systematic effect of about 50 MeV on 
the W mass measurement using the four-jet method, which is of the same order as the total 
measurement error expected in this mode for all four LEP experiments combined. This mass 
shift is found to depend upon the model used and may be larger than estimated above. The 
uncertainty is further increased by other reconnection sources, including perturbative intercon- 
nection, interplay between perturbative and non-perturbative effects [|^, virtual effects, and 



so on. The 50 MeV number above is therefore to be viewed as a realistic uncertainty. 

It is also interesting to study event topologies at LEPl. For instance, a qggq event is nor- 
mally expected to have the partons ordered along a single string, but an alternative subdivision 
into a qq and a gg singlet is possible. Such "reconnection" phenomena (though not quite of 
the same character as the ones we worry about for the W mass issue) could give rapidity gaps 
in events where the quark and antiquark are tagged to lie in the same hemisphere with re- 
spect to the thrust axis Generic event shapes presumably cannot lead to any definite 
conclusions; in part, the information on event shapes is used to tune generators in the first 
place. If a properly tuned generator with "reconnection" included could perform as well as 
the conventional ones, it would provide indications that a given approach is not unreasonable. 
Models including "reconnection" which do not, even with appropriate retuning, describe the 
data are less useful. A recent retuning of Ariadne with "reconnection" by DELPHI shows 
that the quality of the description of LEPl data can even be improved. However, it may be 
entirely fortuitous. 



4.3 Bose-Einstein Effects 



Since the hadronization regions of the W"^ and W overlap, it is natural to assume that some 
coherence effects are present between identical low-momentum bosons stemming from different 
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W's due to Bose-Einstein (BE) correlations. How much such effects would inffuence the W 
mass measurement is difficult to answer. Intuitively, since the BE effect favours production of 
identical bosons close in phase-space, one would expect the softest particles from each W to 
be "dragged" closer to each other. This reduces the momentum of the W's and thus increases 
the measured W mass. However, most of our understanding of the details of multi-particle 
production comes from probabilistic hadronization models where BE correlations are absent. 
Only a few attempts have been made to include such effects and to investigate their influence 
on the W mass measurement. 



One such attempt [B3 used an algorithm (LUBOEI) implemented in the Jetset. This 



algorithm is based on the assumptions that BE effects are local in phase space and do not 
alter the event multiplicity. (Notice that any effects on the multiplicity are already, at least 
partially, accounted for by the tuning of the generators to data.) A re- weighting of events with 
a BE factor is then equivalent, in some approximation, to letting the momenta of the bosons 
produced in the hadronization be shifted somewhat, e.g. to reproduce the difference in two- 
particle correlation functions expected for a source with Gaussian distribution in space-time, 

-^^ = l + Aexp(-g^i?^). (48) 

Here R is the source radius and A is an incoherence parameter. Cbe(<5) and CnoBE(Q) are the 
two-particle correlations as functions of relative four-momenta in a world with and one without 
BE effects, respectively. This is a well-defined generator procedure but, since it is not possible 
to switch off BE in the data, the main challenge of experimental BE studies is to define an 
appropriate no-BE reference sample. With this algorithm, the observed correlations at LEP are 



reproduced using A ~ 1 and R Q.'o im. |68|. One problem with this algorithm is that it does 
not inherently conserve energy and momentum. This has to be done separately by rescaling 
all momenta in an event. Therefore, when this algorithm is used on W'''W~ events, identical 
bosons within each W are shifted closer to each other in phase space resulting in an artificial 
negative shift in the measured W mass even if there is no cross-talk between the W's. In Ref. 



54| this artificial shift was corrected for and an estimate for the true shift in the measured W 
mass was obtained. The result was (AM^^) ^ 100 MeV for 170 GeV center of mass energy. 
This shift increases with energy and decreasing source radius. 



A variation of this algorithm has been investigated where not only the momenta of 
identical particles are shifted closer to each other, but also non-identical particles are shifted 
away from each other, so that the ratio of correlation functions of identical and non-identical 
particles is the same as with the original algorithm. In this way, the energy-momentum non- 
conservation can be made much smaller, minimizing possible dependence on the correction 
procedure. The shift in the W mass is in this case smaller, around 30 MeV. But modifying the 
algorithm slightly, not allowing the shifting of non-identical bosons from different W's, again 
gives a mass shift around 100 MeV. 

In both these cases the shift was calculated on the generator level, assuming that the origin 
of each particle was known. It is clear that an experimentally feasible reconstruction algorithm 
may be more or less sensitive to these BE effects. 
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Another attempt to estimate the influence of BE correlations on the W mass measurement 
was based on a simplified toy model of W decays, in which only decays into two or three 
particles were considered ||6l|. Each of the two W's was treated as a scalar, decaying into 



a heavy particle and one or two pions, according to a relativistic Breit-Wigner. The decay 
amplitudes can thus be written down explicitly, with and without symmetrization with respect 
to exchange of the pions. This allowed the generation of such low-multiplicity events, each with 
two weights corresponding to the symmetrized and non- symmetrized cases, respectively, and 
thus the study of the effect on the reconstructed W mass. The mass of the heavy particle in 
the model was chosen to correspond roughly to the total mass of the hadronic decay products 
of a W boson with one or two pions removed, i.e. 75-79 GeV. 

The model calculations confirm that substantial shifts of the W mass due to the BE sym- 
metrization can occur, and also suggests that the magnitude of the effect can depend signifi- 
cantly on the precise way in which the W mass is reconstructed from the data. In the particular 
case studied, the peak position of the mass distribution was less sensitive to the BE effect than 
the average mass. The interpretation was hampered by severe threshold effects, however, and it 
is in any case not clear how to draw quantitative conclusions valid for the multi-particle decays 
that form the bulk of real hadronic W decays. 

At our present level of understanding, effects of BE correlations on the W-mass measure- 
ment of the order of 100 MeV cannot be excluded, though this number may be viewed as a 
"worst case" estimate. Much more work is needed to develop realistic models and to study the 
sensitivity of different reconstruction algorithms. It may also be possible to find signals, e.g. by 
studying the two-particle correlation function using only pairs where the particles come from 
different reconstructed W's, that could enable us to deduce from data the magnitude of the 
effect. 



4.4 Summary and Outlook 

Several activities are under way, especially on the experimental front, so it would be premature 
to draw any definite conclusions. However, it is not excluded that interconnection and Bose- 
Einstein effects could each contribute a 50 MeV mass shift in the hadronic W'''W~ decay 
channel. Most interconnection models predict a positive mass shift, but we know of no general 
argument why it would have to be so. Furthermore, not all possible interconnection effects 
have been studied so far. In contrast, a positive mass shift from BE effects appears plausible 
on general grounds. 

In spite of the bias towards positive shifts, the conservative approach would be to quote a 
symmetric "theory uncertainty" of the order of ±100 MeV. This is to be viewed as some sort 
of theorist's la error, that is, larger numbers can not be excluded but are less likely. The value 
above refers to studies at 175 GeV. The energy dependence has not yet been well studied, but 
indications are that the uncertainty is not smaller even at the maximum LEP2 energy. We note 
here that the interconnection phenomena are not expected to die out rapidly with increasing 
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energy |Q, since the W resonance is so wide that the two W decays appear "on top of each 
other", in terms of hadronic distance scales, over the full LEP2 energy range. 



Interconnection effects might be found in the data itself, at LEP2 or even by a careful 
analysis of LEPl data. If so, it would be very interesting. However, a non-discovery would 
not prove non-existence: so far, we have failed to find realistic signals for some of the current 
approaches. We could only hope to exclude some extreme scenarios with large values of the 
reconnection probability. Conversely, at our present level of understanding, a discovery would 
not guarantee a unique recipe for how to "correct" the data. This is illustrated by BE effects, 
where it should be fairly straightforward to observe cross-talk at LEP2 in principle, but where 
statistics presumably will not allow a sufficiently precise investigation for the effects on the 
measured W mass to be well defined. 

At the end of this workshop, another model of interconnections appeared [^, with which 



its authors find mass shifts of up to several hundred MeV, i.e. more than observed in previous 
studies. This exemplifies the need for continued theoretical and experimental activity in the 
field to better understand the issues. However, in view of the well-known complexity of hadronic 
final states, it would be unreahstic to expect these studies to give any simple answers. Monte 
Carlo models are all we have so far, and may be the best we can get in the foreseeable future. 
Therefore the main conclusion is clear: one cannot blindly average W mass results in the 
hadronic W'^'W" channel with numbers obtained elsewhere. 



5 Conclusions 

The goal of AMw < 0(50) GeV appears to be achievable at LEP2. The 'worst-case' scenario, 
in which only the qq/i/ channel is used in the direct reconstruction method because of colour 
reconnection in the qqqq channel (see below), gives a total estimated error of AMw = 44 MeV, 
assuming four experiments each collecting 500 pb~^ at y/s = 175 GeV. 

Two methods have been studied in detail. A precise measurement of the e^e~ — > W^W~ 
threshold cross section can be compared to the theoretical prediction as a function of Mw, 
see Fig. |^. The maximum statistical power is obtained by running the collider at an optimal 
energy (y/s)"^*' ~ 2Mw + 0.5 GeV ~ 161 GeV. All three decay channels (qqqq, qq/z^ and lulu) 
can be used, and the total estimated errors are listed in Table ^. Some marginal improvement 
in the systematic error can be expected, but ultimately the measurement is statistics limited. 
An integrated luminosity of at least 50 pb~^ per experiment is required to obtain a precision 
comparable to that expected from the combined Fermilab Tevatron experiments. The total error 
using the threshold cross-section method for four experiments each with 50 pb~^ is estimated 
to be 108 MeV. 

One very attractive feature of the threshold measurement is that it appears to fit in very 
well with the anticipated LEP2 machine schedule, viz. a run of at least 25 pb~^ luminosity at 
^/s = 161 GeV in 1996 followed by the bulk of the luminosity at higher energy. 
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Figure 14: Collider energy dependence of the statistical errors. 



The direct reconstruction method, in wliicli Mw is reconstructed from tlie invariant 
mass of tlie decay products, appears to offer the highest precision. The statistical and sys- 
tematic errors are summarised in Table |l^. More work is needed to investigate different fitting 
procedures, different fragmentation parameters etc., and this may lead to small but important 
reductions in the errors. The major outstanding problem concerns the issues of colour reconnec- 
tion, i.e. non-perturbative strong interactions between the decay products of two hadronically 
decaying W bosons which may significantly distort the invariant mass distributions. This ques- 
tion is currently receiving much attention, but at this time it is impossible to say whether the 
problem can be completely overcome. At best, one may eventually hope to demonstrate that 
the effect on the reconstructed Mw value is either small or under theoretical control, so that 
the results from the qqqq and qqHu channels can be combined. At worst, the precise magnitude 
of the colour reconnection effect may remain unknown, with different models giving different 
mass shifts larger or comparable to the remaining systematic and statistical errors. We have 
therefore chosen to present results both for the qq/z/ channel only (where colour reconnection 
is absent) and for the two channels combined. The estimated total errors are 44 MeV and 
34 MeV respectively, for 500 pb~^ luminosity for each experiment at =175 GeV. 
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The collider energy dependence of the direct reconstruction measurement has also been 



studied. The results are summarised in Fig. 14. There is clearly a 'blind window' between 



161 GeV and 170 GeV where a precision measurement of Mt^ is impossible. In this window the 
threshold cross section method loses sensitivity to Mw, and the direct reconstruction method 
has insufficient statistics. 

Increasing the LEP2 energy from 175 GeV to 192 GeV has an essentially neutral effect on 
the measurement of Mw- Most of the anticipated errors, in both the qqqq and qqlu channels, 
appear to be approximately energy independent .0 Two important energy-dependent systematic 
errors, which are common to all experiments and both decay channels, are from initial state 
radiation and from beam energy calibration. Both are expected to increase slightly with energy, 
although the overall effect is marginal. 
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